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NUCLEAR RECEPTOR LIGANDS AND LIGAND BINDING DOMAINS 

ACKNOWLEDGMENTS 
This invention was suppoaed in part by grams from the National Institutes 
of Health grant number 1 ROl DK43787, and 5 ROl DK 41842. The U.S. 
Government may have rights in ihis invention. 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application claims the benefit of the following provisional applications: 
United States Ser. No. 60/008,540 and 60/008,543. filed December 13, 1995, and 
Ser. No. 60/008,606, filed December 14, 1995. 

INTRODUCTION 

Technical Field 

This invention relates to computational methods for designing ligands that 
bind to nuclear receptors, crystals of nuclear receptors, synthetic ligands of 
nuclear receptors and methods of using synthetic ligands. 

Background 

Nuclear receptors represent a superfamily of proteins that specifically bind 
a physiologically relevant small molecule, such as hormone or vitamin. As a 
result of a molecule binding to a nuclear receptor, the nuclear receptor changes the 
ability of a cell to transcribe DNA, i.e. nuclear receptors modulate the 
transcription of DNA, although they may have transcription independent actions. 
Unlike integral membrane receptors and membrane associated receptors, the 
nuclear receptors reside in either the cytoplasm or nucleus of eukaryotic cells. 
Thus, nuclear receptors comprise a class of intracellular, soluble ligand-regulated 
transcription factors. 

Nuclear receptors include receptors for glucoconicoids (GRs). androgens 
(ARs), mineralocorticoids (MRs), progestins (PRs). estrogens (ERs), thyroid 
hormones (TRs), vitamin D (VDRs), retinoids (RARs and RXRs). The so called 
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"orphan receptors" are also part of the nuclear receptor superfamily. as they are 
structurally homologous to the classic nuclear receptors, such as steroid and 
thyroid receptors. To date, ligands have not been identified with orphan receptors 
but It is likely that small molecule ligands will be discovered in the near fumre for 
this class of transcription factors. Generally, nuclear receptors specifically bmd 
physiologically relevant small molecules with high affinity and apparent Kd's are 
commonly in the 0.01 - 20 nM range, depending on the nuclear receptor/ligand 
pair. 

Development of synthetic ligands that specifically bind to nuclear receptors 
has been largely guided by the trial and error method of drug design despite the 
importance of nuclear receptors in a myriad of physiological processes and 
medical conditions such as hypertension, inflammation, hormone dependent 
cancers (e.g. breast and prostate cancer), modulation of reproductive organ 
modulation, hyperthyroidism, hypercholesterolemia and obesity. Previously, new 
ligands specific for nuclear receptors were discovered in the absence of 
information on the three dimensional strucmre of a nuclear receptor with a bound 
ligand. Before the present invention, researchers were essentially discovering 
nuclear receptor ligands by probing in the dark and without the ability to visualize 
how the amino acids of a nuclear receptor held a ligand in its grasp. 

Consequently, it would be advantageous to devise methods and 
'compositions for reducing the time required to discover ligands to nuclear 
receptors, synthesize such compounds and administer such compounds to 
organisms to modulate physiological processes regulated by nuclear receptors. 

SUMMARY OF THE INVENTION 
The present invention provides for crystals of nuclear receptor ligand 
binding domains with a ligand bound to the ligand binding domain (LBD). The 
crystals of the present invention provide excellent atomic resolution of the amino 
acids that interact with nuclear receptor ligand. especially thyroid receptor ligands. 
The three dimensional model of a nuclear receptor LBD with a ligand bound 
reveals a previously unknown strucmre for nuclear receptors and shows that the 
ligand is bound in a water inaccessible binding cavity of the ligand binding domain 
of the nuclear receptor. 



wo 97/21993 PCT/US96/20778 

The present invention also provides for coniputational methods using three 
dimensional models of nuclear receptors that are based on crystals of nuclear 
receptor LBDs. Generally, the computational method of designmg a nuclear 
receptor ligand determines which amino acid or ammo acids of a nuclear receptor 
LBD interact with a chemical moiety (at least one) of tiie ligand using a three 
dimensional model of a crystallized protein comprising a nuclear receptor LBD 
with a bound ligand, and selecting a chemical modification (at least one) of the 
chemical moiety to produce a second chemical moiety with a structure that either 
decreases or increases an interaction between the interacting amino acid and the 
second chemical moiety compared to the interaction between the interacting amino 
acid and the corresponding chemical moiety on the natural hormone. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a diagram illustrating computational methods for designing 
ligands that interact with nuclear receptors of the nuclear receptor superfamily. 

FIG. 2 is a schematic representation of nuclear receptor structures, 
indicating regions of homology within family members and functions of the 
various domains. 

FIG. 3 shows the aligned amino acid sequences of the ligand binding 
domains of several members of the nuclear receptor superfamily. 

FIG. 4 is a ribbon drawing of the rat TR-ct LBD with secondary structure 
elements labelled. The ligand (magenta) is depicted as a space-filling model. 
Alpha helices and coil conformations are yellow, beta strands are blue. 

FIG. 5 shows two cross-sections of a space-filling model of rat TR-oc 
exposmg the ligand (magenta) tightly packed within the receptor. 

FIG. 6 is a schematic of the ligand binding cavity. Residues which interact 
with the ligand appear approximately at the site of interaction. Hydrogen bonds 
are shown as dashed lines between the bonding partners; distances for each bond 
are listed. Non-bonded contacts are shown as radial spokes which face toward 
interacting atoms. 

FIG. 7 is the distribution of crystallographic lemperamre factors in the 
refined rat TR-a LBD. The distribution is represented as a color gradation 
ranging from less than 15 (dark blue) to greater than 35 (yellow-green). 
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FIG. 8 is a ribbon drawing of the rat TR-a LBD showing the c-terminal 
activation domain to ligand. Residues which comprise the c-terminal activation 
domain (Pro393-Phe405) are depicted as a stick representation. Hydrophobic 
residues, panicularly Phe401 and Phe405 (blue) face mwards toward the ligand. 
Glu403 (red) projects outward into the solvent. 

FIG. 9 is an electrostatic potential surface of the rat TR-a LBD, calculated 
using GRAPH. Negative electrostatic potential is red; positive electrostatic 
potential is blue. The c-terminal activation domain forms a largely hydrophobic 
(white). The Glu403 is presented as a singular patch of negative charge (red). 

FIG. 10 is a diagram comparing agonists and antagonists for several 
nuclear receptors. 

FIG. 11 is the synthetic scheme for preparation of TSl, TS2, TS3, TS4 
and TS5. 

FIG. 12 is the synthetic scheme for preparation of TS6 and TS7. 
FIG. 13 is the synthetic scheme for preparation of TS8. 
FIG. 14 is the synthetic scheme for preparation of TSIO. 
FIG. 15 depicts the chemical structures of several TR ligands. 
FIG. 16 is a graph illustrating competition assays in which T3 and triac 
compete with labeled T3 for binding to human TR-of or human TR-^. 

FIG. 17 depicts a Scatchard analysis of labelled T, binding to TR-a and 

TR-^. 

FIG. 18 is a chart showing the effect of TS-10 on the transcriptional 
regulation of the DR4-ALP reponer gene in the presence or absence of T3 as 
assayed in TRAFal reporter cells. 

FIG. 19 is a chart showing the effect of TS-10 on the transcriptional 
regulation of the DR4-ALP reponer gene in the presence or absence of T3 as 
assayed in TRAF/31 reporter cells. 

FIG. 20 is a chart showing the effect of TS-10 on the transcriptional 
regulation of the DR4-ALP reporter gene in the presence or absence of T3 as 
assayed in HepG2, a liver reporter cell line. 
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FIG. 21 is a panial ribbon drawing of TR-a LBD with T3 in the ligand 
binding cavity. Selected interacting annino acids are labelled, including Iie22L 
Ile222 and Ser260, Ala263, Ile299 and Leu 276. 

FIG. 22 is a partial ribbon drawing of TR-a LBD with T3 and Dimil 
5 superimposed in the ligand binding cavity. Interactions with ne221. Ile222. 

AIa260, He 299 and Leu276 are labelled. 

FIG. 23 is a partial ribbon drawing of TR-a LBD with T3. illustrating the 
three Arginine residues (Arg228, Arg262 and Arg 266 (dark stick figures)) of the 
polar pocket, three water molecules HOH502. HOH503 and HOH504, with 
10 hydrogen bonds indicated by dotted lines. 

FIG. 24 is a partial ribbon drawing of TR-a LBD with tnac, illustrating 
the three Arginine residues (dark stick figures) of the polar pocket, water 
molecules (HOH503. HOH504 and HOH600), with hydrogen bonds indicated by 
dotted lines. 

15 FIG. 25 is a panial ribbon drawing of the TR-a LBD with T3 and iriac 

superimposed in the ligand binding cavity. The drawing shows several interacimg 
amino acid residues in the polar pocket that remain unchanged whether T3 or uiac 
occupies the ligand binding cavity: Arg262, Asnl79, HOH503 and HOH504, and 
Ser277. Both Arg228 and Arg 266 occupy two different positions, depending on 

20 whether T3 or triac is bound. 

FIG. 26 is a stereochemical representation of the TRa LBD with Dimit 

bound. 

APPENDIX 1 is an appendix of references. 

APPENDIX 2 is a chart of amino acids that interact with a TR ligand, for 
25 TR complexed with Dimit, Triac, IpBr2, and T3. 



DETAILED DESCRIPTION OF THE INVENTION 

Introduction 

The present invention provides new methods, particularly computational 
30 methods, and compositions for the generation of nuclear receptor synthetic iigands 

based on the three dimensional strxjcture of nuclear receptors, particularly the 
thyroid receptor (herein referred to as "TR"). Previously, the lack of three 
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dimensional structural information about the ligand binding domain of a nuclear 
receptor thwarted the field of nuclear receptor drug discovery, especially the 
absence of three dimensional structural information relatmg to a nuclear receptor 

with a ligand bound. 

Described herein for the first time are crystals and three dimensional 
strucniral information from a nuclear receptor's ligand binding domain (LBD) with 
a ligand bound. Such crystals offer superior resolution at the atomic level and the 
ability to visualize the coordination of nuclear receptor ligands by amino acids that 
comprise the LBD. The present invention also provides computational methods 
for designing nuclear receptor synthetic ligands using such crystal and three 
dimensional strucmral information to generate synthetic ligands that modulate the 
conformational changes of a nuclear receptor's LBD. Such synthetic ligands can 
be designed using the computational methods described herein and shown, in part, 
in FIG. 1. These computational methods are particularly usehil in designing an 
antagonist or partial agonist to a nuclear receptor, wherein the antagonist or partial 
agonist has an extended moiety that prevents any one of a number of ligand- 
induced molecular events that alter the receptor's influence on the regulation of 
gene expression, such as preventing the normal coordination of the activation 
domain observed for a naturally occurring ligand or other ligands that mimic the 
naturally occurring ligand, such as an agonist. As described herein, synthetic 
ligands of nuclear receptors will be useful in modulating nuclear receptor activity 
in a variety of medical conditions. 

Applicability to Nuclear receptors 

The present invention, paaicularly the computational methods, can be used 
to design drugs for a variety of nuclear receptors, such as receptors for 
glucocorticoids (GRs), androgens (ARs), mineralocorticoids (MRs), progestins 
(PRs), estrogens (ERs), thyroid hormones (TRs), vitamin D (VDRs), retinoid 
(RARs and RXRs) and peroxisomal proliferators (PPAP). The present invention 
can also be applied to the "orphan receptors," as they are structurally homologous 
in terms of modular domains and primary structure to classic nuclear receptors, 
such as steroid and thyroid receptors. The ammo acid homologies of orphan 
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receptors with other nuclear receptors ranges from very low (< 15%) to in the 
range of 35% when compared to rat RARa and human TR-/? receptors, for 
example. In addition, as is revealed by the X-ray crystallographic structure of the 
TR and structural analysis disclosed herein, the overall folding of liganded 
superfamily members is likely to be similar. Although ligands have not been 
identified with orphan receptors, once such ligands are identified one skilled in the 
art will be able to apply the present invention to the design and use of such 
ligands, as their overall structural modular motif will be similar to other nuclear 
receptors described herein. 

Modular Functional Domains Of Nuclear receptors 

The present invention will usually be applicable to all nuclear receptors, as 
discussed herein, in part, to the patterns of nuclear receptor activation, structure 
and modulation that have emerged as a consequence of determining the three 
dimensional structures of nuclear receptors with different ligands bound, notably 
the three dimensional structures or crystallized protein structure of the ligand 
binding domains for TR-a and TR-|3. Proteins of the nuclear receptor superfamily 
display substantial regions of amino acid homology, as described herein and 
known in the an see FIG. 2. Members of this family display an overall strucmral 
motif of three modular domains (which is similar to the TR three modular domain 
motiO: 

1) a variable amino-terminal domain; 

2) a highly conserved DNA-binding domain (DBD); and 

3) a less conserved carboxyl-ierminal ligand-binding domain (LBD). 
The modularity of this superfamily permits different domains of each protein to 
separately accomplish different functions, although the domains can influence each 
other. The separate function of a domain is usually preserved when a particular 
domain is isolated from the remainder of the protein. Using conventional protein 
chemistry techniques a modular domain can sometimes be separated from the 
parent protein. Using conventional molecular biology techniques each domain can 
usually be separately expressed with its original function intact or chimerics of two 
different nuclear receptors can be constructed, wherein the chimerics retain the 
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properties of the individual functional domains of the respective nuclear receptors 
from which the chimerics were generated. 

FIG. 2 provides a schematic representation of family member strucmres. 
indicating regions of homology within family members and functions of the 
various domains. 



Amino Terminal Domain 

The amino terminal domain is the least conserved of the three domains and 
varies markedly m size among nuclear receptor superfamily members. For 

10 example, this domain contains 24 amino acids in the VDR and 603 amino acids in 

the MR. This domain is involved in transcriptional activation and in some cases 
its uniqueness may dictate selective receptor-DNA binding and activation of target 
genes by specific receptor isoforms. This domain can display synergistic and 
antagonistic interactions with the domains of the LBD. For example, studies with 

15 mutated and/or deleted receptors show positive cooperativity of the amino and 

carboxy terminal domains. In some cases, deletion of either of these domains will 
abolish the receptor's transcriptional activation functions. 



DN A- Binding Domain 

20 The DBD is the most conserved structure in the nuclear receptor 

superfamily. It usually contains about 70 amino acids that fold into two zinc 
finger motifs, wherein a zinc ion coordinates four cysteines. DBDs contain two 
perpendicularly oriented a-helixes that extend from the base of the first and second 
zinc fingers. The two zinc fingers function in concert along with non-zinc finger 

25 residues to direct nuclear receptors to specific target sites on DNA and to align 

receptor homodimer or heterodimer interfaces. Various amino acids in DBD 
influence spacing between two half-sites (usually comprised of six nucleotides) for 
receptor dimer binding. For example, GR subfamily and ER homodimers bind to 
half-sites spaced by three nucleotides and oriented as palindromes. The optimal 

30 spacings facilitate cooperative interactions between DBDs, and D box residues are 

part of the dimerization interface. Other regions of the DBD facilitate DNA- 
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protein and proiein-protein interactions required for RXR homodimerization and 
heterodimerization on direct repeal elements. 

The LBD may influence the DNA bindmg of the DBD. and the innuence 
can also be regulated by ligand binding. For example, TR ligand binding 
inHuences the degree to which a TR binds to DNA as a monomer or dimer. Such 
dimenzation also depends on the spacmg and orientation of the DNA half sites 

The nuclear receptor superfamily has been subdivided into two subfamilies: 
1) GR (GR, AR, MR and PR) and 2) TR (TR, VDR, RAR, RXR, and most 
orphan receptors) on the basis of DBD stmctures. interactions with heat shock 
proteins (hsp), and ability to form heterodimers. GR subgroup members are 
tightly bound by hsp in the absence of ligand. dimenze following ligand binding 
and dissociation of hsp, and show homology m the DNA half sites to which they 
bind. These half sites also tend to be arranged as palindromes. TR subgroup 
members tend to be bound to DNA or other chromatin molecules when 
unliganded, can bind to DNA as monomers and dimers, but tend to form 
heterodimers, and bind DNA elements with a variety of orientations and spacings 
of the half sites, and also show homology with respect to the nucleotide sequences 
of the half sites. ER does not belong to either subfamily, since it resembles the 
GR subfamily in hsp interactions, and the TR subfamily in nuclear localization and 
DNA-binding propenies. 

Ligand Binding Domain 

The LBD is the second most highly conserved domain in these receptors. 
Whereas integrity of several different LBD sub-domains is important for ligand 
binding, truncated molecules containing only the LBD retain normal ligand-binding 
activity. This domain also panicipates in other functions, including dimerizaiion. 
nuclear translocation and transcriptional activation, as described herein. 
Importantly, this domain binds the ligand and undergoes ligand-induced 
conformational changes as detailed herein. 

Most members of the superfamily. including orphan receptors, possess at 
least two transcription activation subdomains, one of which is constimiive and 
resides in the amino terminal domain (AF-1). and the other of which (AF-2 (also 
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referenced as TAU 4)) resides in the ligand-binding domain whose activity is 
regulated by binding of an agonist ligand. The function of AF-2 requires an 
activation domain (also called transactivation domain) that is highly conserved 
among the receptor superfamily (approximately amino acids 1005 to 1022). Most 
LBDs contain an activation domain. Some mutations in this domam abolish AF-2 
function, but leave ligand binding and other functions unaffected. Ligand binding 
allows the activation domain to serve as an interaction site for essential co- 
activator proteins that function to stimulate (or in some cases, inhibit) 
transcription. 

The carboxy-terminal activation subdomain, as described herein is in close 
three dimensional proximity in the LBD to the ligand, so as to allow for ligands 
bound to the LBD to coordinate (or mteract) with amino acid(s) in the activation 
subdomain. As described herein, the LBD of a nuclear receptor can be expressed, 
crystallized, its three dimensional structure determined with a ligand bound (either 
using crystal data from the same receptor or a different receptor or a combination 
thereoO, and computational methods used to design ligands to its LBD, 
particularly ligands that contain an extension moiety that coordinates the activation 
domain of the nuclear receptor. 

Once a computationally designed ligand (CDL) is synthesized as described 
herein and known in the art, it can be tested using assays to establish its activity as 
an agonist, partial agonist or antagonist, and affinity, as described herein. After 
such testing, the CDLs can be further refmed by generating LBD crystals with a 
CDL bound to the LBD. The strucmre of the CDL can then be further refmed 
using the chemical modification methods described herein for three dimensional 
models to improve the activity or affmity of the CDL and make second generation 
CDLs with improved properties, such as that of a super agonist or antagonist 
described herein. 

Nuclear Receptor Isoforms 

The present invention also is applicable to generatmg new synthetic ligands 
to distinguish nuclear receptor isoforms. As described herein, CDLs can be 
generated that distinguish between isoforms. thereby allowing the generation of 
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either tissue specific or function specific synthetic ligands. For instance, GR 
subfamily members have usually one receptor encoded by a single gene, with the 
exception that there are two PR isoforms, A and B, translated from the same 
mRNA by alternate initiation from different AUG codons. This method is 
especially applicable to the TR subfamily which usually has several receptors that 
are encoded by two (TR) or three (RAR, RXR, and PPAR) genes or have alternate 
RNA splicing and such an example for TR is described herein. 



Nuclear Receptor Crystals 

10 The invention provides for crystals made from nuclear receptor ligand 

binding domains with the ligand bound to the receptor. As exemplified in the 
Examples, TRs are crystallized with a ligand bound to it. Crystals are made from 
purified nuclear receptor LBDs that are usually expressed by a cell culture, 
preferably E. coli. Preferably, different crystals (co-crystals) for the same nuclear 

15 receptor are separately made using different ligands, such as a naturally occurring 

ligand and at least one bromo- or iodo- substituted synthetic ligand that acts as an 
analog or antagonist of the naturally occurring ligand. Such bromo- and iodo- 
substitutions act as heavy atom substitutions in nuclear receptor ligands and 
crystals of nuclear receptor proteins. This method has the advantage for phasing 

20 of the crystal in that it bypasses the need for obtaining traditional heavy metal 

derivatives. After the three dimensional structure is determined for the nuclear 
receptor LED with its ligand bound, the three dimensional structure can be used in 
computational methods to design a synthetic ligand for the nuclear receptor and 
further activity structure relationships can be determined through routine testing 

25 using the assays described herein and known in the an. 

Expression and Purification of other Nuclear Receptor LED Structures 

High level expression of nuclear receptor LBDs can be obtained by the 
techniques described herein as well as others described in the literature. High 
30 level expression in E. coli of ligand binding domains of TR and other nuclear 

receptors, including members of the steroid/thyroid receptor superfamily, such as 
the estrogen (ER). androgen (AR). mineralocorticoid (MR), progesterone (PR). 
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RAR, RXR and vitamin D (VDR) receptors can also be achieved. Yeast and other 
eukaryotic expression systems can be used with nuclear receptors that bind heat 
shock proteins as these nuclear receptors are generally more difficult to express in 
bacteria, with the exception of ER, which can be expressed in bacteria. 
Representative nuclear receptors or their ligand binding domains have been cloned 
and sequenced; human RAR-a, human RAR-7, human RXR-a, human RXR-/3, 
human PPAR-a, human PPAR-^, human PPAR-7, human VDR, human ER (as 
described in Seielstad et al.. Molecular Endocrinology, vol 9:647-658 (1995), 
incorporated herein by reference), human GR, human PR, human MR, and human 
AR. The ligand binding domain of each of these nuclear receptors has been 
identified and is shown in FIG. 3. Using the information in FIG. 3 in 
conjunction with the methods described herein and known in the art, one of 
ordinary skill in the an could express and purify LBDs of any of the nuclear 
receptors, including those illustrated in FIG. 3, bind it to an appropriate ligand, 
and ciystailize the nuclear receptor's LBD with a bound ligand. 

FIG. 3 is an alignment of several members of the steroid/thyroid hormone 
receptor superfamily that indicates the amino acids 10 be included in a suitable 
expression vector. 

Extracts of expressing cells are a suitable source of receptor for purification 
and preparation of crystals of the chosen receptor. To obtain such expression, a 
vector is constructed in a manner similar to that employed for expression of the rat 
TR alpha (Apriletti et aL Protein Expression and Purification, 6:368-370 (1995), 
herein incorporated by reference). The nucleotides encoding the amino acids 
encompassing the ligand binding domain of the receptor to be expressed, for 
example the estrogen receptor ligand binding domain (hER-LBD) (corresponding 
to R at position 725 to L at position 1025 as standardly aligned as shown in the 
FIG. 3), are inserted into an expression vector such as the one employed by 
Apriletti et al (1995). For the purposes of obtaining material that will yield good 
crystals it is preferable to include at least the amino acids corresponding to human 
TR-/3 positions 725 to 1025. Stretches of adjacent amino acid sequences may be 
included if more structural information is desired. Thus, an expression vector for 
the human estrogen receptor can be made by insening nucleotides encoding amino 
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acids from position 700 to the c-terminus at position 1071. Such a vector gives 
high yield of receptor in E. coli that can bind hormone (Seielstad et al. Molecular 
Endocrinology Vol 9:647-658 (1995)). However, the c-terminal region beyond 
position 1025 is subject to variable proteolysis and can advantageously be excluded 
from the constmct, this technique of avoiding variable proteolysis can also be 
applied to other nuclear receptors. 

TR-ct And TR-^ As Examples of Nuclear receptor LBD Structure and Function 
TR Expression, Purification And Crystallization 

As an example of nuclear receptor structure of the ligand binding domain 
the a- and ^3- isoforms of TR are crystallized from proteins expressed from 
expression constructs, preferably constructs that can be expressed in E. coli. 
Other expression systems, such as yeast or other eukaryotic expression systems 
can be used. For the TR, the LBD can be expressed without any portion of the 
DBD or amino-terminal domain. Portions of the DBD or amino-terminus can be 
included if further strucairal information with amino acids adjacent the LBD is 
desired. Generally, for the TR the LBD used for crystals will be less than 300 
amino acids in length. Preferably, the TR LBD will be at least 150 amino acids in 
length, more preferably at least 200 amino acids in length, and most preferably at 
least 250 amino acids in length. For example the LBD used for crystallization can 
comprise amino acids spanning from Met 122 to Val 410 of the rat TR-a, Glu 202 
to Asp 461 of the human TR-^. 

Typically TR LBDs are purified to homogeneity for crystallization. Purity 
of TR LBDs is measured with SDS-PAGE, mass spectrometry and hydrophobic 
HPLC. The purified TR for crystallization should be at least 97.5 % pure or 
97.5%, preferably at least 99.0% pure or 99.0% pure, more preferably at least 
99.5% pure or 99.5% pure. 

Initially purification of the unliganded receptor can be obtained by 
conventional techniques, such as hydrophobic interaction chromatography (HPLC). 
ion exchange chromatography (HPLC). and heparin affinity chromatography. 

To achieve higher purification for improved crystals of nuclear receptors, 
especially the TR subfamily and TR. it will be desirable to ligand shift purify the 
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nuclear receptor using a column that separates the receptor according to charge, 
such as an ion exchange or hydrophobic interaction column, and then bind the 
eluted receptor with a hgand, especially an agonist. The ligand induces a change 
in the receptor's surface charge such that when re-chromatographed on the same 
column, the receptor then elutes at the position of the liganded receptor are 
removed by the original column run with the unliganded receptor, Usually 
saturating concentrations of ligand are used in the column and the protein can be 
preincubated with the ligand prior to passing it over the column. The structural 
studies detailed herein indicate the general applicability of this technique for 
obtaining super-pure nuclear receptor LBDs for crystallization. 

More recently developed methods involve engineering a "tag" such as with 
histidine placed on the end of the protein, such as on the amino terminus, and then 
using a nickJe chelation column for purification, Janknecht R., Proc. Natl Acad. 
Sci. USA Vol 88:8972-8976 (1991) incorporated by reference. 

To determine the three dimensional structure of a TR LBD, or a LBD from 
another member of the nuclear receptor superfamily, it is desirable to co-crystalize 
the LBD with a corresponding LBD ligand. In the case of TR LBD, it is 
preferable to separately co-crystaiize it with ligands such as T3, IpBr and Dimit 
that differ in the heavy atoms which they contain. Other TR ligands such as those 
encompassed by Formula 1 described herem and known in the prior an. can also 
"be used for the generation of co-crystals of TR LBD and TR ligands. Of the 
compounds encompassed by Formula 1 it is generally desirable to use at least one 
ligand that has at least one bromo- or iodo- substitution at the R,, R5, Ry or R5' 
position, preferably such compounds will be have at least two such substitutions 
and more preferably at least 3 such substitutions. As described herein, such 
substitutions are advantageously used as heavy atoms to help solve the phase 
problem for the three dimensional structure of the TR LBD and can be used as a 
generalized method of phasing usmg a halogen (e.g. I or Br) substituted ligand, 
especially for nuclear receptors. 

Typically purified LBD, such as TR LBD, is equilibrated at a saturating 
concentration of ligand at a temperature that preserves the integrity of the protein. 
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Lieand equilibration can be established between 2 and ST'' C. although the 
receptor tends to be more stable in the 2-20*' C range. 

Preferably crystals are made with the hanging drop methods detailed 
herein. Regulated temperature control is desirable to improve crystal stability and 
quality. Temperatures between 4 and 25°C are generally used and it is often 
preferable to test crystallization over a range of temperatures. In the case of TR it 
is preferable to use crystallization temperatures from 18 to 25°C, more preferably 
20 to 23^*0, and most preferably 22°C. 

Complexes of the TR-a LBD with a variety of agonists, including 3,5,3'- 
triiodothyronine (T3), 3*-isopropyl-3,5-dibromothyronine (IpBr.), 3'-isopropyl-3.5- 
dimethylthyronine (Dimit), and 3,5,3'-triiodothyroacetic acid (triac), are prepared 
with by methods described herein. Cocrystals of the rTR-a LBD, with ligand 
prebound, are prepared by vapor diffusion at ambient temperature from 15% 2- 
methyl-2,4-penlanediol (MPD). The crystals are radiation sensitive, and require 
freezing to measure complete diffraction data. On a rotating anode X-ray source, 
the crystals diffract to --3A; synchrotron radiation extends the resolution limit 
significantly, to as high as 2.0A for T3 cocrystals. The composition of the thyroid 
hormone, combined with the ability to prepare and cocrystallize the receptor 
complexed with a variety of analogs, permitted the unusual phasing strategy. This 
phasing strategy can be applied to the ligands of the nuclear receptors described 
therein by generating I and Br substitutions of such ligands. In this strategy, 
cocrystals of the TR LBD containing four hormone analogs that differ at the 3,5, 
and 3* positions (T3, IpBrj, Dimit, and triac) provided isomorphous derivatives. 
For this set of analogs, the halogen substituents (2Br and 31 atoms) function as 
heavy atoms, while the Dimit cocrystal (3 alkyl groups) acts as the parent. The 
initial 2.5A multiple isomorphous replacement/anomalous scattering/density 
modified electron density map allowed the LBD to be traced from skeletons 
created in the molecular graphics program 05 (Jones, T.A. et al. ACTA Cryst, 
47:110-119 (1991), incorporated by reference herein). A model of the LBD was 
built in four fragments, Argl57-Glyl84, Trpl86-Gly 197, Serl99-Pro205, and 
VaI210-Phe405, and refined in XPLOR using positional refinement and simulated 
annealing protocols. Missing residues were built with the aid of difference 
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density. The final model was refined to R^^ys. = 21.8% and R^,, = 24.4% for 
data fromJ5.0 to 2.2A, see Table 3. 

This phasing strategy can be applied to the ligands of the nuclear receptors 
described herein by generating I and Br substirutions of such Ugands. 

Three Dimensional Structure of TR LBD 
Architecture of TR LBD 

As an example of the three dimensional structure of a nuclear receptor, the 
folding of the TR-a, LBD is shown in FIG. 4. The TR-a LBD consists of a single 
structural domain packed in three layers, composed of twelve a-helices, Hl-12, 
and four shon /3-strands, Sl-4, forming a mixed jS-sheet. The buried hormone and 
three amiparallel a-helices, H5-6, H9, and HIO, form the central layer of the 
domain, as shown in FIG. 4. HI, H2. H3 and SI form one face of ihe LBD, 
with the opposite face formed by H7, H8, HIL and H12. The first 35 amino 
acids of the N-terminus (Metl22-Glnl56) are not visible in the electron density 
maps. The three dimensional strucmre of the heterodimeric RXR:TR DNA- 
binding domains bound to DNA, amino acids Met 122 - Glnl51 of the TR DBD 
make extensive contacts with the minor groove of the DNA8. The five disordered 
amino acids (Argl52-Glnl56), which reside between the last visible residue of the 
TR DBD and the first visible residue of the LBD likely represent the effective 
"hinge" linking the LBD and the DBD in the intact receptor. 

The predominantly helical composition and the layered arrangement of 
secondary structure is identical to that of the unliganded hRXRa, confirming the 
existence of a common nuclear receptor fold between two nuclear receptors. 

The TR LBD is visible beginning at Argl57, and continues in an extended 
coil conformation to the start of HI. A turn of a-helix, H2, covers the hormone 
binding cavity, immediately followed by short /3-strand, SI, which forms the edge 
of the mixed /3-sheet, parallel to S4. the outermost of the three antiparallel strands. 
The chain is mostly irregular until H3 begins, amiparallel to HI. H3 bends at 
Ile221 and IIe222, residues which contact the ligand. The chain turns almost 90' 
at the end of H3 to form an incomplete a-helix. H4. The first buried core helix, 
H5^6, follows, its axis altered by a kink near the ligand at Gly 253. The helix is 
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composed of mostly hydrophobic sidechams mtemipted by iwo sinking exceptions: 
Arg262 is solvent inaccessible and interacts with the ligand carboxylate (1- 
subsntuent), and Glu256 meets Arg329 from H9 and Arg375 from Hll m a polar 
invagination. H5-6 termmates in a short i3-sirand. S2. of the four strand mixed 
sheet. S3 and S4 are joined through a left-handed turn, and further linked by a 
salt bridge between Lys284 and Asp272. Following S4. H7 and H8 form an L, 
stabilized by a salt bridge between Lys268 and Asp277. The mm between H7 and 
H8 adopts an unusual conformation, a result of interaction with ligand and us 
glycine rich sequence. H9 is the second core helix, antiparailel to the neighboring 
H5-6. Again, two buried polar sidechains are found, Glu315 and Gln320. 
Glu315 forms a buried salt bridge with His358 and Arg356. The oxygen of 
Gln320 forms a hydrogen bond with the buried sidechain of His 175. The chain 
then switches back again to form HIO, also antiparailel to H9. HU extends 
diagonally across the full length of the molecule. Immediately after Hll, the chain 
forms a type II mm, at approximately 90* to Hll. The chain then mms again to 
form H 12, which packs loosely against H3 and Hll as part of the hormone or 
ligand binding cavity. The final five amino acids at the C^terminus, Glu406 - 
Val410, are disordered. 

TR LBD's Ligand Binding Cavity As An Example Of A Nuclear Receptor's 

Buried Ligand Cavity 

The three dimensional sirucmre of the TR LBD leads to the startling 
finding that ligand binding cavity of the LBD is solvent inaccessible when a T3 or 
its isostere is bound to the LBD. This surprising result leads to a new model of 
nuclear receptor three dimensional strucmre and function, as further described 
herein, particularly in the sections elucidating the computational methods of ligand 
design and the application of such methods to designing nuclear receptor synthetic 
ligands that contain extended positions that prevent normal activation of the 
activation domain. 

Dimit, the ligand bound to the receptor, is an isostere of T^ and a thyroid 
hormone agonist. Therefore the binding of Dimit should refiect that of T„ and the 
Dimit-bound receptor is expected to be the active conformation of TR. The ligand 



wo 97/21993 



PCT/L'S96/20778 



18 



10 



15 



20 



25 



is buried within the receptor, providing the hydrophobic core for a subdomain of 
the protein, as shown in FIG. 5 a and b. H5-6 and H9 comprise the hydrophobic 
core for the rest of the receptor. 

An extensive binding cavity is constructed from several stnjctural elements. 
The cavity is enclosed from above by H5-6 (Met 256- Arg266), from below by H7 
and H8 and the intervening loop (Leu287- Ile299), and along the sides by H2 
(185-187), by the turn between S3 and S4 (Leu276-Ser277), by H3 (Phe215- 
Arg228), by HU (His381-Met388) and by H12 (Phe401-Phe405). The volume of 
the cavity defined by these elements, calculated by GRASP (Columbia University, 
USA) (600 A3), is essentially the volume of the hormone (530 A). The remaining 
volume is occupied by water molecules surrounding the amino-propionic acid 
substituent. FIG. 6 depicts various contacts (or interactions) between TR's LBD 
and the ligand. 

The planes of the inner and outer (prime ring) rings of the ligand are 
rotated from planarity about 60° with respect to each other, adopting the 3'-disial 
conformation (in which the 3' substituent of the outer ring projects down and away 
from the inner ring). The ammo-propionic acid and the outer phenolic ring 
assume the transoid conformation, each on opposite sides of the irmer ring. The 
torsion angle Xi for the amino- propionic acid is 300 . 

The amino-propionic acid substituent is packed loosely in a polar pocket 
formed by side chains from H2, H4 and S3. The carboxylate group forms direct 
hydrogen bonds with the guanidium group of Arg228 and the amino N of Ser277. 
In addition. Arg262, Arg266 and Asnl79 interact with the carboxylate through 
water-mediated hydrogen bonds. The three arginine residues create a significantly 
positive local electrostatic potential, which may stabilize the negative charge of the 
carboxylate. No hydrogen bond is formed by the amino nitrogen. The 
interactions of the amino-propionic acid substiment are consistent with the fact that 
triac. which lacks the amino nitrogen, has a binding affmity equal to that of Tj, 
indicating that the amino nitrogen and longer aliphatic chain of T3 do not 
contribute greatly to binding affinity. 

The diphenyl ether, in contrast, is found buried within the hydrophobic 
core. The inner ring packs in a hydrophobic pocket formed by H3, H5-6, and S3. 
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Pockets for the 3- and 5-methyl substituents are not completely filled, as expected 
since the van der waals radius of methyl substituent for Dimit is smaller than the 
lodme substituent provided by the thyroid hormone T3. Such pockets are typically 
25 to 100 cubic angstroms (although smaller pocket for substitutes are 
contemplated in the 40 to 80 cubic angstrom range) and could be filled more 
lightly with better fitting chemical substitutions, as described herein. 

The outer ring packed tightly in a pocket formed by H3, H5-6, H7, H8, 
H 11 and HI 2, and the loop between H7 and H8. The ether oxygen is found in a 
hydrophobic environment defined by Phe218, Leu287. Leu276. and L£u292. The 
absence of a hydrogen bond to the ether oxygen is consistent with its role in 
establishing the correct stereochemistry of the phenyl rings, as suggested by potent 
binding of hormone analogs with structurally similar linkages possessing reduced 
or negligible hydrogen bonding capability. The 3'-isopropyl substituent contacts 
Gly290 and 291. The presence of glycine at this position in the pocket can 
explain the observed relationship between activity and the size of 3'-substituents. 
Activity is highest for 3'-isopropyU and decreases with added bulk. The only 
hydrogen bond in the hydrophobic cavity is formed between the phenolic hydroxyl 
and His381 Ne2. The conformation of His381 is stabilized by packing contacts 
provided by Phe405, and Met256. 

The presence of a 5' substituent larger than hydrogen affects the binding 
affinity for hormone. The more abundant thyroid hormone, 3,5,3\5*-tetraiodo-L- 
thyronine (T4), contains an iodine at this position, and binds the receptor with 2% 
of the affinity of T,. The structure suggests that discrimination against is 
accomplished through the combination of steric conflict by Met256 and possibly 
the constraints imposed by the geometry of the hydrogen bond from His381 to the 
phenolic hydroxyl. The 5' position is a preferred location for introducing a 
chemical modification of C-H at the 5' of T3 or and TR agonist, as described 
herein, that produces an extension from the prime ring and results in the creation 
of an antagonist or partial agonist. 

Deletion and antibody competition smdies suggest the involvement of 
residues Prol62 to Val202 in ligand binding. The region does not directly contact 
hormone in the bound stnicmre, although H2 packs against residues forming the 
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polar pocket that interacts with the amino-propionic acid group. One role for H2, 
then, is to stabilize these residues in the bound state, H2, with /3-strands S3 and 
S4, might also represent a prevalent entry point for ligand, since the amino- 
propionic acid of the ligand is oriented toward this region. Studies of receptor 
binding to T3 affinity matrices demonstrate that only a linkage to the amino- 
propionic acid is tolerated, suggesting that steric hindrance present in other 
hnkages prevent binding. Furthermore, the crystallographic temperature factors 
suggest the coil and /J-strand region is most flexible part of the domain FIG. 7. 
Participation of this region, part of the hinge domain between the DBD and LBD, 
in binding hormone may provide structural means for ligand binding to influence 
DNA binding, since parts of the Hinge domain contact DNA. 

TR LBD Transcriptional Activation Helix As An Example Of A Nuclear Receptor 
Activation Domain 

In addition to the startling fmding that the ligand binding cavity is solvent 
inaccessible when loaded with a ligand, the activation helix of TR LBD presents a 
surface to the ligand cavity for interaction between at least one amino acid and the 
bound ligand. The C-ierminal 17 amino acids of the TR, referred to as the 
activation helix or AF-2 (an example of an LBD activation domain), are implicated 
in mediating hormone-dependent transcriptional activation. Although, mutations of 
key residues within the domain decrease ligand-dependent activation it was unclear 
until the present invention whether such mutations directly affected ligand 
coordination. Although some mutations of this domain have been noted to reduce 
or abolish ligand binding, other mutations in more distant sites of the LBD have a 
similar effect. 

Activation domains among nuclear receptors display an analogous three 
dimensional relationship to the binding cavity, which is a region of the LBD that 
binds the molecular recognition domain of a ligand, i.e. the activation domain 
presents a portion of itself to the binding cavity (but necessarily the molecular 
recognition domain of the ligand). Many nuclear receptors are expected to have 
such domains, including the retinoid receptors. RAR and RXR, the glucocorticoid 
receptor GR, and the estrogen receptor ER. Based upon the TR's sequence, the 
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domain is proposed to adopt an amphipathic helical strucaire. ^-sheet or mixed 
secondary structures, could be present as activaiion domains in less related nuclear 
receptors. 

Within the activation domain, the highly conserved motif 4'<I'XE4><{>, where 
4> represents a hydrophobic residue, is proposed to mediate interactions between 
the receptors and transcriptional coactivators. Several proteins have been 
identified which bind the TR in a hormone-dependent fashion. One of these. 
Tripl, is related to a putative yeast coactivator SugU and also interacts with both 
the C-terminal activation domain and a subset of the basal transcriptional 
machinery, suggesting a role in transactivation by the TR. Other proteins, such as 
R1P140, SRCl, (Onate, S.A. et. aK. Science 270:1354-1357 (1995)) and TF-1 (see 
also Ledouarim, B., et. al., EMBO J. 14:2020-2033 (1995)), also interact with 
other nuclear receptors in a Ugand dependent manner through the C-terminal 
domain. Binding of these proteins can be modulated using the TR ligands 
described herein especially those TR ligands with extensions that stericaily hinder 
the interaction between the highly conserved motif and other proteins. 

The C-terminal activation domain of the TR forms an amphipathic helix, 
H12, which nestles loosely against the receptor to form part of the hormone 
binding cavity. The helix packs with the hydrophobic residues facing inward 
towards the hormone binding cavity, and the charged residues, including the 
highly-conserved glutamate, extending into the solvent, as shown in FIG. 8. The 
activation helix of TR LBD presents Phe 401 to the ligand binding cavity and 
permits direct coordination with the hormone i.e. such amino acids interact with 
the ligand forming a van der waals contact with the plane of the outer phenyl ring. 
Phe 405 also interacts with His 381, perhaps stabilizing its hydrogen bonding 
conformation, i.e. a favorable hydrogen bond interaction. Participation of Phe 
401 and Phe 405 in binding hormone explains how mutation of these residues 
decreases hormone binding affinity. Furthermore, the impact of these mutations 
on activation likely derives from a role in stabilizing the domain in the bound 
structure through increased hydrogen bond interaction of dipole interactions. Giu 
403 extends into the solvent, emphasizing its critical role in transactivation. In its 
observed conformation, presented on the surface as an ordered residue, against a 
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background of predominantly hydrophobic surface, Glu 403 is available to interact 
with activator proteins described herein, as shown in FIG. 9. The other charged 
residues, Glu 405 and Asp 406 are disordered, as the helix frays at Phe 405, 

Two other sequences in the TR, t2 and r3. activate transcription when 
expressed as fiision proteins with a DNA-binding domain. The sequences, 
discovered in the TRB, correspond to TR-a residues Prol58-Ilel68 in HI (t2), 
and Gly290-Leu3 19 in H8 and H9 (t3). Unlike the C-terminal activation domain. 
r2 and t3 do not appear to represent modular structural units in the rat TR-Of 
LBD, nor present a surface for protein-protein interactions: the critical 
aspanate/glutamate residues of r3 are located on two separate helices, and do not 
form a single surface; the charged residues of r2 are engaged in ion pair 
interactions with residues of the LBD. Thus, t2 and t3 may not function as 
activation domains in the context of the entire receptor. 

Computational Methods For Designing A Nuclear Receptor LBD LiGAND 

The elucidation of the three dimensional strucoire of a nuclear receptor 
ligand binding domain provides an important and useful approach for designing 
ligands to nuclear receptors using the computational methods described herein. By 
inspecting the FIGLFRES it can be determined that the nuclear receptor ligand is 
bound in a water inaccessible binding cavity in the LBD and that chemical 
moieties can be added to selected positions on the ligand. Such chemical 
modifications, usually extensions, can fill up the binding cavity represented in the 
FIGURES for a tighter fit (or less water) or can be used to disrupt or make 
contacts with amino acids not in contact with the ligand before the chemical 
modification was introduced or represented in a figure of the three dimensional 
model of the LBD. Ligands that interact with nuclear superfamily members can 
act as agonists, antagonists and partial agonists based on what ligand^nduced 
conformational changes take place. 

Agonists induce changes in receptors that place them in an active 
conformation that allows them to influence transcription, either positively or 
negatively. There may be several different ligand-induced changes in the 
receptor's conformation. 
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Antagonists, bind to receptors, but fail to induce conformational changes 
that alter the receptor's transcriptional regulatory properties or physiologically 
relevant conformations. Binding of an antagonist can also block the binding and 
therefore the actions of an agonist. 

Partial agonists bind to receptors and induce only part of the changes in the 
receptors that are induced by agonists. The differences can be qualitative or 
quantitative. Thus, a partial agonist may induce some of the conformation 
changes induced by agonists, but not others, or it may only induce certain changes 
to a limited extent, 

Ligand-induced Conformational Changes 

As described herein, the unliganded receptor is in a configuration that is 
either inactive, has some activity or has repressor activity. Binding of agonist 
ligands induces conformational changes in the receptor such that the receptor 
becomes more active, either to stimulate or repress the expression of genes. The 
receptors may also have non-genomic actions, some of the known types of 
changes and/or the sequelae of these are listed herein. 
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Heat Shock Protein Binding 

For many of the nuclear recepiors ligand binding induces a dissociation of 
heat shock proteins such that the receptors can form dimers in most cases, after 
which the receptors bind to DNA and regulate transcription. 

Nuclear receptors usually have heat shock protein binding domains that 
present a region for binding to the LBD and can be modulated by the binding of a 
ligand to the LBD. Consequently, an extended chemical moiety (or more) from 
the ligand that stabilizes the binding or contact of the heat shock protein binding 
domain with the LBD can be designed using the computational methods described 
herein to produce a partial agonist or antagonist. Typically such extended 
chemical moieties will extend past and away from the molecular recognition 
domain on the ligand and usually past the buried binding cavity of the ligand. 

Dimerization and Heterodimerization 

With the receptors that are associated with the hsp in the absence of the 
ligand, dissociation of the hsp results in dimerization of the receptors. 
Dimerization is due to receptor domains in both the DBD and the LBD. Although 
the main stimulus for dimerization is dissociation of the hsp. the ligand- induced 
conformational changes in the receptors may have an additional facilitative 
influence. With the receptors that are not associated with hsp in the absence of 
the ligand, particularly with the TR, ligand binding can affect the pattern of 
dimerization/heterodimerization. The influence depends on the DNA binding site 
context, and may also depend on the promoter context with respect to other 
proteins that may interact with the receptors. A common pattern is to discourage 
monomer formation, with a resulting preference for heterodimer formation over 
dimer formation on DNA. 

Nuclear receptor LBDs usually have dimerization domams that present a 
region for binding to another nuclear receptor and can be modulated by the 
binding of a ligand to the LBD. Consequently, an extended chemical moiety (or 
more) from the ligand that disrupts the binding or contact of the dimerization 
domain can be designed using the computational methods described herein to 
produce a partial agonist or antagonist. Typically such extended chemical moieties 
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will extend past and away from the molecular recogmtion domam on the ligand 
and usually past the buried binding cavity of the ligand. 

DNA Binding 

In nuclear receptors that bind lo hsp, the ligand-induced dissociation of hsp 
with consequent dimer formation allows, and therefore, promotes DNA bindmg. 
With receptors that are not associated (as in the absence of ligand), ligand binding 
tends to stimulate DNA binding of heterodimers and dimers, and to discourage 
monomer binding to DNA. However, with DNA contaimng only a single half 
site, the ligand tends to stimulate the receptor's binding to DNA. The effects are 
modest and depend on the nature of the DNA site and probably on the presence of 
other proteins that may interact with the receptors. Nuclear receptors usually have 
DBD (DNA binding domains) that present a region for binding to DNA and this 
binding can be modulated by the binding of a ligand to the LBD. Consequently, 
an extended chemical moiety (or more) from the ligand that disnipts the binding or 
contact of the DBD can be designed using the computational methods described 
herein to produce a partial agonist or antagonist. Typically such extended 
chemical moieties will extend past and away from the molecular recognition 
domain on the ligand and usually past the buried binding cavity of the ligand. 

Repressor Binding 

Receptors that are not associated with hsp in the absence of ligand 
frequently act as transcriptional repressors in the absence of the ligand. This 
appears to be due, in part, to transcriptional repressor proteins that bind to the 
LBD of the receptors. Agonist binding induces a dissociation of these proteins 
from the receptors. This relieves the inhibition of transcription and allows the 
transcriptional transaciivation functions of the receptors to become manifest. 

Transcriptional Transactivation Functions 

Ligand binding induces transcriptional activation functions in two basic 
ways. The first is through dissociation of the hsp from receptors. This 
dissociation, with consequent dimenzation of the receptors and their binding to 
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DNA or other proteins in the nuclear chromatin allows transcriptional regulator)- 
properties of the receptors to be manifest. This may be especially ir\ie of such 
functions on the amino terminus of the receptors. 

The second way is to alter the receptor to mteract with other proteins 
involved in transcription. These could be protems that interact directly or 
indirectly with elements of the proximal promoter or proteins of the proximal 
promoter. Alternatively, the interactions could be through other transcription 
factors that themselves interact directly or indirectly with proteins of the proximal 
promoter. Several different proteins have been described that bind to the receptors 
in a ligand-dependent manner. In addition, it is possible that in some cases, ihc 
ligand-induced conformational changes do not affect the binding of other proteins 
to the receptor, but do affect their abilities to regulate transcription. 

Nuclear receptors or nuclear receptor LBDs usually have activation 
domains that present a region for binding to DNA and can be modulated by the 
binding of a ligand to the LBD. Consequently, an extended chemical moiety (or 
more) from the ligand that disrupts the binding or contact of the activation domain 
can be designed using the computational methods described herein to produce a 
partial agonist or antagonist. Typically such extended chemical moieties will 
extend past and away from the molecular recognition domain on the ligand and 
usually past the buried binding cavity of the ligand and in the direction of the 
activation domain, which is often a helix as seen in the three dimensional model 
shown in the FIGURES in two dimensions on paper or more conveniently on a 
computer screen. 

Ligand-induced Conformational Change 

Plasma proteins bind hormones without undergoing a conformational 
change through a static binding pocket formed between monomers or domains. 
For example, the tetrameric thyroid-binding plasma protein transthyrecm forms a 
solvent-accessible hormone-binding channel at the oligomer interface. The 
strucnire of the protein is unchanged upon binding hormone with respect to the 
appearance of a buried binding cavity with a ligand bound. 
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However, the structural role for a iigand bound to a nuclear receptor LBD, 
like rat TR-a LBD, predicts that the receptor would differ in the bound and 
unbound states. In the absence of hormone, the receptor would possess a cavity at 
its core, uncharacteristic of a globular protein. A ligand (e.g. hormone) completes 
the hydrophobic core of the active receptor after it binds to the nuclear receptor. 
Ligand binding by the receptor is a dynamic process, which regulates receptor 
function by inducing an altered conformation. 

An exact description of the hormone-induced conformational changes 
requires comparison of the structures of the liganded and the unliganded TR. The 
structure of the unliganded human RXRa may substitute as a model for the 
unliganded TR. The rat TR-a LBD and human RXRa LBDs adopt a similar fold, 
and it is likely that the structural similarity extends to the conformational changes 

after ligand binding. 

There are three major differences between the two structures, which indeed 
appear to be the result of ligand binding- First, the bound rat TR-a LBD structure 
is more compact, with the hormone tightly packed within the hydrophobic core of 
the receptor. By contrast, the unliganded human RXRa LBD contains several 
internal hydrophobic cavities. The presence of such cavities is unusual in folded 
proteins, and is likely a reflection of the unJiganded state of the receptor. Two of 
these cavities were proposed as possible binding sites for 9-cis retinoic acid, 
though these multiple sites only panly overlap with the single buried binding 
cavity observed in the liganded rat TR-a LBD. 

The second difference involves Hll in the rat TR-a LBD, which 
contributes part of the hormone binding cavity. Hll. continuous in the rat TR-a 
LBD, is broken at Cys 432 in the RXR, forming a loop between HIO and Hll in 
the hRXRa. This residue corresponds to His381 in the TR, which provides a 
hydrogen bond to the outer ring hydroxy I of the ligand. Furthermore, the 
hormone binding cavity occupied by ligand in the rat TR-a LBD is mterrupted in 
the hRXRa by the same loop, forming an isolated hydrophobic pocket in the RXR 
with H6 and H7. In the bound rat TR-a LBD, the corresponding helices H7 and 
H8 are contiguous with the binding pocket, and enclose the hormone binding 
cavity from below. 



wo 97/21993 



PCT/US96/20778 



28 



10 



15 



20 



25 



The third difference between the two receptors is the position of the C- 
terminal activation domain. While the C-terminal activation domain forms a- 
hehces in both receptors, the domain in the rat TR-a LBD follows a proline-nch 
turn, and lies against the receptor to contribute pan of the binding cavity. In 
contrast, the activation domain in the unliganded hRXRa. is part of a longer helix 
which projects into the solvent. 

These differences lead to a model for an alternate conformation of the TR 
LBD assumed in the absence of ligand. In the unliganded TR, the subdomain of 
the receptor surrounding the hormone binding cavity is loosely packed, with the 
binding cavity occluded by a partly unstructured HU providing a partial core for 
the receptor. 

Upon binding hormone, residues which form a coil in the unbound receptor 
engage the ligand, and continues Hll. The ordering of HIl could unblock the 
hydrophobic cavity, allowing H7 and H8 to interact with hormone. The extended 
hydrophobic cavity then collapses around the hormone, generating the compact 
bound structure. 

It is possible to predict iigand-induced conformational changes in the C- 
terminal activation domain that rely, in part, on an extended structure in the 
unliganded TR that repacks upon ligand binding. The ligand- induced 
conformation change can be subtle since the amino acid sequence of the rat TR-a 
in the mm (393-PTELFPP-399) significantly reduces the propensity of the peptide 
chain of the rat TR-a to form an a-helix and therefore repacking can be 
accomplished with a minor change in volume. 

After the ligand-induced conformational change occurs, it is likely that the 
conformation of the C-terminal activation domain in the bound strucmre changes 
packing compared to the unbound form of the receptor. Binding of the ligand 
improves the stability of the activation domain. The activation domain packs 
loosely even in the bound structure, as measured by the distribution of packing 
interactions for the entire LBD. The packing density for the activation domain, 
defined as the number of atoms within 4.5A, is 1.5 standard deviations below the 
mean. For comparison, another surface helix, HI, is 0.5 standard deviations 
below the mean and the most poorly packed pan of the stnicture. the irregular coil 



PCT/TJS96/20778 

WO 97/21993 

29. 

from residues Ilel96-Asp206, is 2.0 standard deviations below the mean. 
Moreover, the majority of packing contacts for the C-terminal domain m the 
bound receptor are provided either by residues which interact with ligand. such as 
His38K or by the ligand itself. The conformation of these residues can be 
expected to be different in the bound and unbound receptors, and by extension the 
conformation of C-terminal activation domain which relies upon these interactions. 
Without the stabilization provided by a bound ligand, it is likely that the C- 
terminal activation domain is disordered prior to hormone binding. 

The interrelation of ligand- induced conformational changes is evident as 
described herein. For example, His381 from Hll and Phe405 from H12 interact 
in the bound structure to provide a specific hydrogen bond to the phenolic 
hydroxyl. The ligand-induced changes which affect Hll and Hi2 are remforcing. 
and lead to the formation of the compact, bound state. 

Computational Methods Using Three Dimensional Models and 
Extensions Of Ligands 

The three-dimensional structure of the liganded TR receptor is 
unprecedented, and will greatly aid in the development of new nuclear receptor 
synthetic ligands, such as thyroid receptor antagonists. In addition, this receptor 
superfamily is overall well suited to modem methods including three-dimensional 
structure elucidation and combinatorial chemistry such as those disclosed in EP 
335 628, U.S. patent 5,463,564, which are incorporated herein by reference. 
Structure determination using X-ray crystallography is possible because of the 
solubility properties of the receptors. Computer programs that use crystallography 
data when practicing the present invention will enable the rational design of ligand 
to these receptors. Programs such as RASMOL can be used with the atomic 
coordinates from crystals generated by practicing the invention or used to practice 
the invention by generating three dimensional models and/or determining the 
structures involved in ligand binding. Computer programs such as INSIGHT and 
GRASP allow for further manipulation and the ability to introduce new strucmres. 
In addition, high throughput binding and bioactivity assays can be devised using 
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purified recombinant protein and modem reponer gene transcription assays 
described herein and known in the an in order to refme the activity of a CDL. 

Generally the computational method of designing a nuclear receptor 
synthetic ligand comprises two steps: 

1) determining which amino acid or amino acids of a nuclear receptor LBD 
interacts with a first chemical moiety (at least one) of the ligand using a three 
dimensional model of a crystallized protein comprising a nuclear receptor LBD 
with a bound ligand, and 

2) selecting a chemical modification (at least one) of the first chemical 
moiety to produce a second chemical moiety with a structure to either decrease or 
increase an interaction between the interacting amino acid and the second chemical 
moiety compared to the interaction between the mteracting amino acid and the first 
chemical moiety. 

As shown herein, interacting amino acids form contacts with the ligand and 
the center of the atoms of the interacting amino acids are usually 2 to 4 angstroms 
away from the center of the atoms of the ligand. Generally these distances are 
determined by computer as discussed herein and in McRee 1993, however 
distances can be determined manually once the three dimensional model is made. 
Examples of interacting amino acids are described in Appendix 2. See also 
Wagner et ah. Nature 378(6558);670-697 (1995) for stereochemical figures of 
three dimensional models. More commonly, the atoms of the ligand and the 
atoms of interacting amino acids are 3 to 4 angstroms apart. The invention can be 
practiced by repeating steps 1 and 2 to refine the fit of the ligand to the LBD and 
to determine a better ligand, such as an agonist. As shown in the FIGURES the 
three dimensional model of TR can be represented in two dimensions to determine 
which amino acids contact the ligand and to select a position on the ligand for 
chemical modification and changing the interaction with a particular amino acid 
compared to that before chemical modification. The chemical modification may 
be made using a computer, manually using a two dimensional representation of the 
three dimensional model or by chemically synthesizing the ligand. The three 
dimensional model may be made using Appendix 2 and the FIGURES. As an 
additional step, the three dimensional model may be made using atomic 
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coordinates of nuclear receptor LBDs from crystallized protein as known tn the 




art, see McRee 1993 referenced herein. 




The ligand can also interact with distant ammo acids after chemical 




modification of the ligand to create a new ligand. Distant ammo acids are 
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generally not in contact with the ligand before chemical modification. A chemical 
modification can change the structure of the ligand to make as new ligand that 
interacts with a distant amino acid usually at least 4.5 angstroms away from the 
ligand. Often distant amino acids will not line the surface of the binding cavity 
for the ligand, as they are too far away from the ligand to be pan of a pocket or 


10 


surface of the binding cavity. 

The interaction between an atom of a LBD amino acid and an atom of an 
LBD ligand can be made by any force or attraction described in nature. Usually 
the interaction between the atom of the amino acid and the hgand will be the result 
of a hydrogen bonding interaction, charge interaction, hydrophobic interaction. 


15 


van der waals interaction or dipole interaction. In the case of the hydrophobic 
interaction it is recognized that this is not a per se interaction between the amino 
acid and ligand, but rather the usual result, in part, of the repulsion of water or 
other hydrophilic group from a hydrophobic surface. Reducing or enhancing the 
interaction of the LBD and a ligand can be measured by calculating or testing 


20 


binding energies, computationally or using thermodynamic or kinetic methods as 
known in the art. 

Chemical modifications will often enhance or reduce interactions of an 
atom of a LBD amino acid and an atom of an LBD ligand. Steric hmderancc will 
be a common means of changing the interaction of the LBD binding cavity with 


25 


the activation domain. Chemical modifications are preferably introduced at C-H, 
C- and C-OH position in ligands, where the carbon is part of the ligand structure 
which remains the same after modification is complete. In the case of C-H, C 
could have 1, 2 or 3 hydrogens, but usually only one hydrogen will be replaced. 
The H or OH are removed after modification is complete and replaced with the 


30 


desired chemical moiety. 

Because the thyroid receptor is a member of the larger superfamily of 
hormone-binding nuclear receptors, the rules for agonist and antagonist 
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development will be recognized by one skilled in the art as useful in designing 
ligands to the entire superfamily. Examining the structures of ioiown agonists and 
antagonists of the estrogen and androgen receptors supports the generality of 
antagonist mechanism of action as shown in FIG. 10. 

The overall folding of the receptor based on a comparison of the reported 
structure of the unliganded RXR and with amino acid sequences of other 
superfamily members reveals that the overall folding of receptors of the 
superfamily is similar. Thus, it is predicted from the structure that there is a 
general pattern of folding of the nuclear receptor around the agonist or antagonist 
ligand. 

The three dimensional structure of a nuclear receptor with a ligand bound 
leads to the nonobvious observation that a nuclear receptor folds around agonist 
ligands, as the binding cavity fits the agonist, especially the agonist's molecular 
recognition domain, and antagonists commonly have chemical structures that 
extend beyond the ligand, especially the agonist, and would prohibit folding of the 
receptor around the ligand to form a buried binding cavity or other groups that 
have the same effect. The location of the extension could affect the folding in 
various ways as indicated by the strucmre. Such extensions on antagonists are 
shown in FIG. 10 for various receptors and compared to the corresponding 
agonist. 

For example, an extension towards the carboxy-terminal activation helix 
affects the packing/folding of this helix into the body of the receptor. This in turn 
can affect the ability of this ponion of the nuclear receptor to mteract with other 
proteins or other portions of the receptor, including transcriptional transactivalion 
functions on the opposite end of the linear receptor, or the receptor's amino 
terminus that may interact directly or indirectly with the carboxy-terminal 
transactivation domain (including helix 12). Extensions in this direction can also 
affect the packing of helix 11 of TR (or its analogous helix in nuclear receptors) 
into the body of the receptor and selectively affect dimerization and 
heterodimerization of receptors. An extension pointing towards helix 1 can affect 
the relationship of the DNA binding domain and hinge regions of the receptors 
with the ligand binding domain and selectively or in addition affect the receptors' 
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binding to DNA and/or interactions of receptors with proteins that interact with 
this region of the receptor. Other extensions towards helix 11 can be made to 
affect the packing of this helix and helices 1 and 10 and thereby dimerizaiion. 
Such chemical modifications can be assessed using the computational methods 
described herein. It is also possible that, in some cases, extensions may protrude 
through the receptor that is otherwise completely or incompletely folded around 
the ligand. Such protruding extensions could present a steric blockade to 
mteractions with co-activators or other proteins. 

The three dimensional strucmre with the ligand buried in the binding cavity 
immediately offers a simple description of a nuclear receptor that has a binding 
cavity that contains hinges and a lid, composed of one or more structural 
elements, that move to accommodate and surround the ligand. The ligand to TR 
can be modified on specific sites with specific classes of chemical groups that will 
serve to leave the lid and hinge region in open, partially open or closed states to 
achieve partial agonist or antagonist functions. In these states, the biological 
response of the TR is different and so the structure can be used to design 
particular compounds with desired effects. 

Knowledge of the three-dimensional structure of the TR-T3 complex leads 
to a general model for agonist and antagonist design. An important novel feauire 
of the structural data is the fact that the T, ligand is completely buried within the 
central hydrophobic core of the protein. Other ligand-receptor complexes 
belonging to the nuclear receptor superfamily will have a similarly buried ligand 
binding site and therefore this model will be useful for agonistyantagonist design 
for the entire superfamily. 

When design of an antagonist is desired, one needs either to preserve the 
important binding contacts of namral hormone agonist while incorporating an 
"extension group" that interferes with the normal operation of the ligand-receptor 
complex or to generate the requisite binding affinity through the interactions of the 
extensions with receptor domains. 

The model applied to antagonist design and described herein is called the 
"Extension Model." Antagonist compounds for nuclear receptors should contain 
the same or similar groups that facilitate high-affinity binding to the receptor, and 
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in addition, such compounds should contain a side chain which may be large 
and/or polar. This side chain could be an actual extension, giving it bulk, or it 
could be a side group with a charge function that differs from the agonist iigand 
For example, substitution of a CH, for CH.OH at the 21 -position, and alteration at 
the 1 1 -position from an OH group to a keto group of conisol generates 
glucocorticoid antagonist activity (Robsseau, G,G., et. ai., J. Mol. Biol. 67:99- 
115 (1972)). However, in most cases effective antagonists have more bulky 
extensions. Thus, the antiglucocorticoid (and antiprogestin) RU486 contains a 
bulky side group at the 11-position (Horwitz, K.B. Endocrine Rev. 13: 146-163 
(1992)). The antagonist compound will then bind within the buried iigand binding 
site of the receptor with reasonably high affmity (100 nM), but the extension 
function will prevent the receptor- Iigand complex from adopting the necessary 
conformation needed for transcription factor function. The antagonism (which 
could be in an agonist or antagonist) may manifest itself at the molecular level in a 
number of ways, including by preventing receptor homo/heterodimer formation at 
the HRE, by preventing coactivator binding to receptor monomers, homodimers or 
homo/heterodimers, or by a combination of these effects which otherwise prevent 
transcription of hormone responsive genes mediated by ligand-induced effects on 
the HRE. There are several antagonist compounds for nuclear receptors in the 
prior art (see also Horwitz, K.B., Endocrine Rev. 13:146-163 (1992), Raunnaud 
J.P. et. al., J. Steroid Biochem. 25:811-833 (1986). Keiel S., et. al.. Mol. Cell. 
Biol. 14:287-298 (1994) whose antagonist function can be explained by the 
extension hypothesis. These compounds are shown in FIG. 10 along with their 
agonist counterparts. Each of these antagonists contains a large extension group 
attached to an agonist or agonist analogue core structure. Importantly, these 
antagonist compounds were discovered by chance and not designed with a 
structure -function hypothesis such as the extension principle. 

One method of design of a thyroid antagonist using the extension 
hypothesis is provided below as a teaching example. The three-dimensional 
structure of the TR-aDimit complex combined with structure-activity data 
published in the prior art, especially those reference herein, can be used to 
establish the following iigand-receptor interactions which are most critical for 
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high-affinity ligand binding. A physical picaire of these interactions is shown jn 
FIG. 6. The figure describes the isolated essential contacts for ligand binding. 
Because the ligand is buried in the center of the receptor, the sirucmral spacing 
between these isolated interactions is also important. Thus, our present knowledge 
of this system dictates that, for this example, a newly designed ligand for the 
receptor must contain a thyronine structural skeleton, or two subsiinited aryl 
groups joined by a one-atom spacer. 

The general structure for an antagonist designed by the extension 
hypothesis is exemplified m the following general description of the substituents of 
a TR antagonist (referring to Formula 1): Rl can have anionic groups such as a 
carboxylaie, phosphonate, phosphate, sulfate or sulfite and is connected to the ring 
with a 0 to 3 atom linker, comprising one or more C, O, N, S atoms, and 
preferably a 2 carbon linker. Such Rl can be optionally substituted with an amine 
(e.g. -NH2). R3 and R5 are small hydrophobic groups such as -Br, -I, or -CH3. 
R3 and R5 can be the same substiments or different. R^' can be a hydrophobic 
group that may be larger than those of R3 and R5, such as -I, -CH3, -isopropyl, - 
phenyl, -benzyl, 5 and 6 ring heterocycles. R,* is a group that can participate in a 
hydrogen bond as either a donor or acceptor. Such groups include -OH, -NHj. 
and -SH. R5' is an important extension group that makes this compound an 
antagonist. R5' can be a long chain alkyl (e.g. 1 to 9 carbons, straight chain or 
branched), aiyl (benzyl, phenyl and substimted benzyl and phenyl rings (e.g with 
halogen, alkyl (1 and 5 carbons) and optionally connected to the ring by a -CH2-). 
hctcrocycle (e.g. 5 or 6 atoms, preferably 5 carbons and 1 nitrogen, or five 
carbons), which can optionally include polar (e.g. -OH, -NHj. and -SH), cationic 
(e.g. -NH3, N(CH)3), or anionic (carboxylate, phosphonate. phosphate or sulfate) 
groups. Rs^ can also be a polar (e.g. -OH. -NH^, and -SH), cationic (e.g. -NH3, - 
N(CH3)3), and anionic (carboxylate, phosphonate, phosphate or sulfate) groups. 
X is the spacer group that appropriately positions the two aromatic rings. This 
group is usually a one-atom spacer, such as O, S, SO, S02, NH, NZ where Z is 
an alkyl, CH2, CHOH, CO, C(CH3)0H, and C(CH3)(CH3). R2. R6, R2' and 
R6' can be -F, and -CI and are preferably H. 
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A TR ligand can also be described as a substiaited phenylated 3,5 diiodo 
tyrosine with substimted R5' and R3' groups. R5' can be a long chain alkyl (e.g. 
4 to 9 carbons, straight chain or branched), aryl (benzyl, phenyl and substituted 
benzyl and phenyl rings (e.g with halogen, alkyl (1 and 5 carbons) and optionally 
connected to the ring by a -CH2-), heterocycle (e.g. 5 or 6 atoms, preferably 5 
carbons and 1 nitrogen, or five carbons), which can optionally include polar (e.g. 
-OH. -NHj, and -SH), cationic (e.g. -NH3, N(CH)3), or anionic (carboxylate, 
phosphonate, phosphate or sulfate) groups. R5' can also be a polar (e.g. -OH, 
-NH,, and -SH), cationic (e.g. -NH3, N(CH)3), and aniomc (carboxylate. 
phosphonate, phosphate or sulfate) groups. R3' can be -IsoPr, halogen, -CH3. 
alkyl (1 to 6 carbons) or aryl (benzyl, phenyl and substimted benzyl and phenyl 
rings (e.g with halogen, alkyl (1 and 5 carbons) and optionally connected to the 
ring by a -CH2-), heterocycle (e.g. 5 or 6 atoms, preferably 5 carbons and 1 
nitrogen, or five carbons), which can optionally include polar (e.g. -OH, -NH^. 
and -SH), cationic (e.g. -NH3, N(CH)3), or anionic (carboxylate, phosphonate, 
phosphate or sulfate) groups. 

A TR antagonist can also be a modified agonist (having a diphenyl 
structure) wherein R5' is alkyl, aryl, 5- or 6-membered heterocyclic aromatic, 
heteroalkyi, heteroaryl, arylalkyl, heteroaryl alkyl, polyaromatic, 
polyheteroaromatic, polar or charged groups, wherein said R5' may be substituted 
with polar or charged groups. The R5' groups are defined, as described herein. 

Using these methods the ligands of this example preferably have the 

following properties: 

1. The compounds should bind to the TR with high affmity (for example 

100 nM). 

2. The compounds should bind the receptor in the same basic orientation 
as the naairal hormone. 

3. The extension group R5' should project toward the activation helix 
(C-terminal helix) of the receptor, 

4. The appropriate substiment at R5' should perturb the activation heiix 
from its optima! local structure needed for mediatmg transcription. 
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Antagonists may also be designed with multiple extensions in order to 
block more than one aspect of the folding at any time. 

TR ligands (e.g. super agonists) can be designed (and synthesized) to 
enhance the mteraction of at least one amino acid with at least one chemical 
moiety on the iigand's molecular recognition domain. One method is to enhance 
the charge and polar interactions by replacing the carboxylate of (Rl position) 
with phosphonate, phosphate, sulfate or sulfite. This enhances the interaction with 
Arg 262, Arg 266 and Arg 228. The interaction of at least one amino acid with at 
least one chemical moiety on the Iigand's molecular recognition domain can also 
be enhanced by increasing the size of Rl group to fill the space occupied by water 
when Dimit is bound (referring to Rl). Preferably the group has a 
complementary charge and hydrophobicity to the binding cavity. 

Another way of improving the interaction of at least one amino acid with at 
least one chemical moiety on the iigand's molecular recognition domain is to 
restrict the conformation of the dihedral angle between the two phenyl rings of the 
thyronine hgand in solution. In solution the planes of two phenyl rings are 
orthogonal where the dihedral angle is 90°. In the TR Dimit structure, the 
dihedral angle is close to 60^ A TR ligand design that fixes the angle between 
the two phenyl rings will lead to tighter binding. Such a ligand may be made by 
connecting the R6* and the R5 positions of a thyronine or a substituted thyronme- 
like diphenyl. The size of the cyclic connection can fix the angle between the two 
phenyl rings. Referring specifically to Formula 1, the following cyclic 
modifications are preferred: 1) R, is connected to 2) R3 is connected to R,' or 
3) R5 is connected to R«' and R3 is connected to R2\ The connections can be 
made by an alky! or heteroalkyl chain having between 1 to 6 atoms and preferably 
from 2 to 4 carbon atoms or other atoms. Any position of the heteroalkyl chain 
can be N, O, P or S. The S and P heteroaioms along said heteroalkyl chain are in 
any of their possible oxidative slates. The N heteroatom or any carbon along the 
alkyl or heteroalkyl chain may have one or more Z substiments, wherein Z is 
alkyl, heteroalkyl, aryl, heteroaryl, 5- or 6-membered heterocyclic aromatic. 
These compounds can be claimed with the proviso that Formula 1 does not include 
any prior art compound as of the filing date of this application. 
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The interaction of at least one amino acid with at least one chemical moiety 
on the ligand's molecular recognition domain can also be enhanced by selecting a 
chemical modification that fills the unfilled space between a TR ligand and the 
LBD in the area of the bridging oxygen (such as in T3, iriac or Dimit). Thus, a 
slighter larger moiety that replaces the ether oxygen can enhance binding. Such a 
linker may be a mono- or geminal- disubstituted carbon group A group 
approximately the same size as oxygen but with greater hydrophobicity is 
preferred as well as small, hydrophobic groups for the disubstituted carbon. 

TR-ff AND TR-i3 Selectivity for the Thyroid Hormone Receptor 

Using the method described herein ligands can be designed that selectively 
bind 10 the alpha more than the beta TR. The X-ray crystallographic structure of 
the rat TR-a LBD provides insight into design of such ligands. 

The three dimensional structure reveals that the major difference between 
the TR-a and TR-/3 in the ligand binding cavity resides in amino acid Ser 277 
(with the side group -CH20H) in the rat TR-a and whose corresponding residue is 
331, asparagine (with the side group -CH2CONH2), in the human TR-i3. The side 
chain in human TR-^ is larger, charged and has a different hydrogen bonding 
potential, which would allow the synthesis of compounds that discriminate between 
this difference. 

For example, in the complex of TRa with triac, Ser277 does not participate 
in ligand binding. The absence of a role for Ser277 (Asn331 in beta) is consistent 
with the equal affinity of triac for the alpha and beta isoforms, and indirectly 
suppons the contention that alpha/beta selectivity resides in the amino acid 
substimtion Ser277 to Asn331 and its interaction with Arg228. 

In terms of ligand design, these differences mean that for /3-selective 
ligands, some or all of the following differences should be exploited: 

1. The presence of a larger side chain asparagine. 

2. The ability of the carbonyl group on the side chain to provide a 
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strong hydrogen bond acceptor. 

The ability of the amido group on the side chain to provide a two 
hydrogen bond donors. 
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4. Adjustment of polarity to reorganize the trapped water in the T3 
pocket. 

In terms of pharmaceutical design, these differences mean that for 
a-selective ligands. some or all of the following differences should be exploited: 
I. The presence of a smaller side group. 

2 The ability of the hydroxy 1 on the -CH20H side group carbonyl 
group on the side chain to provide a weak hydrogen donor. 

3. Adjustment of polarity to reorganize the trapped water in the T3 
pocket. 

In both cases these differences can be exploited in a number of ways. For 
example, they can also be used with a software set for construction of novel 
organic molecules such as LUDI from Biosym-MSI. 



Methods of Treatment 
15 The compounds of Formula I can be useful in medical treatments and 

exhibit biological activity which can be demonstrated in the following tests: 

(i) the induction of mitochondrial a-glycerophosphate dehydrogenase 

(GPDH:EC 1.1.99.5). This assay is particularly useful since in certain species 

e.g. rats it is induced specifically by thyroid hormones and thyromimetics in a 
20 close-related manner in responsive tissues e.g. liver, kidney and the heart 

(Westerfield, W.W., Richert, D.A. and Ruegamer, W.R., Endocrinology, 1965. 

77, 802). The assay allows direct measurement in rates of a thyroid hormone-like 

effect of compounds and in particular allows measurement of the direct thyroid 

hormone-like effect on the hean; 
25 (ii) the elevation of basal metabolic rate as measured by the increase in 

whole body oxygen consumption (see e.g.. Barker et al., Ann. N. Y. Acad. Sci., 

86:545-562 (I960)); 

(iii) the stimulation of the rate of beating of atria isolated from animals 
previously dosed with thyromimetrics (see e.g., Stephan et al., Biochem. 
30 Pharmacol- (1992) 13:1969-1974; Yokoyama et aL. J. Med. Chem. 38:695-707 

(1995)); 
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(iv) the change in total plasma cholesterol levels as determined usmg a 
cholesterol oxidase kit (for example, the Merck CHOD iodine colorimetric kit. 
see also, Stephan et al. (1992)); 

(v) the measurement of LDL (low density lipoprotein) and HDL (high 
density lipoprotein) cholesterol in lipoprotein fractions separated by 
ultracentrifugation; and p (vi) the change in total plasma triglyceride levels as 
determined using enzymatic color tests, for example the Merck System GPO-PAP 
method. 

The compounds of Formula 1 can be found to exhibit selective 
thyromimetic activity in these tests, 

(a) by increasing the metabolic rate of test animals, and raising hepatic 
GPDH levels at doses which do not significantly modify cardiac GPDH levels. 

(b) by lowering plasma cholesterol and triglyceride levels, and the ratio 
of LDL to HDL cholesterol at doses which do not significantly modify cardiac 
GPDH levels. 

The compounds of Formula 1 may therefore be used in therapy, in the 
treatment of conditions which can be alleviated by compounds which selectively 
mimic the effects of thyroid hormones in certain tissues whilst having little or no 
direct thyromimetic effect on the heart. For example, compounds of Formula 1 
which raise hepatic GPDH levels and metabolic rate at doses which do not 
significantly modify cardiac GPDH levels are indicated in the treatment of obesity. 

Agonists of Formula 1 will lower total plasma cholesterol, the ratio of 
LDL-cholesterol to HDL-cholesterol and triglyceride levels at doses which do not 
significantly modify cardiac GPDH levels are indicated for use as general 
antihyperlipidaemic (antihyperlipoproteinaemic) agents i.e. in the treatment of 
patients having elevated plasma lipid (cholesterol and triglyceride) levels. In 
addition, in view of this effect on plasma cholesterol and triglyceride, they are 
also indicated for use as specific antt-hypercholesterolemic and anti- 
hypertriglyceridaemic agents. 

Patients having elevated plasma lipid levels are considered at risk of 
developing coronary heart disease or other manifestations of atherosclerosis as a 
result of their high plasma cholesterol and/or triglyceride concentrations. Further, 
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since LDL^cholestero! is believed to be the lipoprotein which induces 
atherosclerosis, and HDL-cholesterol believed to transport cholesterol from blood 
vessel walls to the liver and to prevent the build up of atherosclerotic plaque, anti- 
hyperlipidemic agents which lower the ratio of LDL-cholesterol to HDL 
cholesterol are indicated as anti-atherosclerotic agents, herein incorporated by 
reference U.S. patents 4,826,876 and 5.466.861. 

The present invention also provides a method of producing selective 
thyromimetic activity in certain tissues except the heart which comprises 
administering to an animal in need thereof an effective amount to produce said 
activity of a compound of Formula 1 or a pharmaceutically acceptable salt thereof. 

The present invention also relates to a method of lowering plasma lipid 
levels and a method of lowering the ratio of LDL-cholesierol to HDL-cholesterol 
levels by suitably administering a compound of this invention or a 
pharmaceutically acceptable sale thereof. 

In addition, compounds of Formula 1 may be indicated in thyroid hormone 
replacement therapy in patients with compromised cardiac function. 

In therapeutic use the compounds of the present invention are usually 
administered in a standard pharmaceutical composition. 

The present invention therefore provides in a further aspect pharmaceutical 
compositions comprising a compound of Formula 1 or a pharmaceutically 
acceptable salt thereof and a pharmaceutically acceptable carrier. Such 
compositions include those suitable for oral, parenteral or rectal administration. 

Pharmaceutical Compositions 

Compounds of Formula 1 and their pharmaceutically acceptable salts which 
are active when given orally can be formulated as liquids for example syrups, 
suspensions or emulsions, tablets, capsules and lozenges. 

A liquid composition will generally consist of a suspension or solution of 
the compound or pharmaceutically acceptable salt in a suitable liquid carrier(s). for 
example ethanol, glycerine, sorbitol, non-aqueous solvent such as polyethylene 
glycol, oils or water, with a suspending agent, preservative, surfactant, wetting 
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agent, flavoring or coloring agent. Alternatively, a liquid formulation can be 
prepared from a reconstimtable powder. 

For example a powder containing active compound, suspending agent, 
sucrose and a sweetener can be reconstituted with water to form a suspension; and 
a syrup can be prepared from a powder containing active ingredient, sucrose and a 
sweetener. 

A composition in the form of a tablet can be prepared using any suitable 
pharmaceutical carrier(s) routinely used for preparing solid compositions. 
Examples of such carriers include magnesium stearate, starch, lactose, sucrose, 
microcrystalline cellulose and binders, for example polyvinylpyrrolidone. The 
tablet can also be provided with a color film coating, or color included as part of 
the carrier(s). In addition, active compound can be formulated in a controlled 
release dosage form as a tablet comprising a hydrophilic or hydrophobic matrix. 

A composition in the form of a capsule can be prepared using routine 
encapsulation procedures, for example by incorporation of active compound and 
excipients into a hard gelatin capsule. Alternatively, a semi-solid matrix of active 
compound and high molecular weight polyethylene glycol can be prepared and 
filled into a hard gelatin capsule; or a solution of active compound in polyethylene 
glycol or a suspension in edijble oil, for example liquid paraffm or fractionated 
coconut oil can be prepared and filled into a soft gelatin capsule. Compound of 
Formula 1 and their pharmaceutical ly acceptable salts which are active when given 
parenterally can be formulated for intramuscular or intravenous administration. 

A typical composition for intra-muscular administration will consist of a 
suspension or solution of active ingredient in an oil. for example arachis oil or 
sesame oil. A typical composition for intravenous administration will consist of a 
sterile isotonic aqueous solution containing, for example active ingredient, 
dextrose, sodium chloride, a co-solvent, for example polyethylene glycol and, 
optionally, a chelating agent, for example ethylenediamine tetracetlc acid and an 
anti-oxidant, for example, sodium metabisulphiie. Alternatively, the solution can 
be freeze dried and then reconstituted with a suitable solvent just prior to 
administration. 
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Compounds of stnicture (1) and their phannaceancally acceptable salts 
which are active on rectal administration can be formulated as suppositories. A 
typical suppository formulation will generally consist of active ingredient with a 
binding and/or lubricating agent such as a gelatin or cocoa butter or other low 
melting vegetable or synthetic wax or fat. 

Compounds of Formula 1 and their pharmaceuticaily acceptable salts which 
are active on topical administration can be formulated as transdermal 
compositions. Such compositions include, for example, a backing, active 
compound reservoir, a control membrane, liner and contact adhesive. 

The typical daily dose of a compound of Formula 1 varies according to 
individual needs, the condition to be treated and with the route of administration. 
Suitable doses are in the general range of from 0.001 to 10 mg/kg bodyweight of 

the recipient per day. 

Within this general dosage range, doses can be chosen at which the 
compounds of Formula 1 lower plasma cholesterol levels and raise metabolic rate 
with little or no direct effect on the hean. In general, but not exclusively, such 
doses will be in the range of from 0.5 to 10 mg/kg. 

In addition, within the general dose range, doses can be chosen at which 
the compounds of Formula I lower plasma cholesterol levels and have little or no 
effect on the heart without raising metabolic rate. In general, but not exclusively, 
such doses will be in the range of from 0.001 to 0.5 mg/kg. 

It is to be understood that the 2 sub ranges noted above are not mutually 
exclusive and that the particular activity encountered at a particular dose will 
depend on the nature of the compound of Formula I used. 

Preferably, the compound of Formula 1 is in unit dosage form, for 
example, a tablet or a capsule so that the patient may self-administer a single dose. 
In general, umt doses contain in the range of from 0.05-100 mg of a compound of 
Fonnula 1. Prefen-ed umt doses contain from 0.05 to 10 mg of a compound of 
Formula 1. 

The active ingredient may be administered from 1 to 6 times a day. Thus 
daily doses are in general in the range of from 0.05 to 600 mg per day. 
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Preferably, daily doses are m the range of from 0.05 to 100 mg per day. Most 
preferably from 0.05 to 5 mg per day. 



EXAMPLES 

5 Example 1 - Synthesis of TR Ligands 

Many TR ligands are kjiown in the art, including T4 (thyroxine), T3, T2 
and TS-9. See Jorgensen, Thyroid Hormones and Analogs, in 6 Hormonal 
Proteins and Peptides, Thyroid Hormones 107-204 (Choh Hao Li ed.. 1978), 
incorporated by reference herein.. 

10 The syntheses of several TR ligands are described below. 

Synthesis of TSl, TS2, TS3, TS4, TS5 

TSi, TS2, TS3, TS4 and TS5 and analogs thereof can all be prepared by 
simple acylation of the nitrogen atom of any thyronine analog, includmg T3 

15 (3,5,3'-tri!odo-L-thyronine), T4 (thyroxine) and 3,5-diiodothyronine. TSl and 

TS2 are synthesized by reacting T3 with PhjCHCO.NHS (N-hydroxy succinimide- 
2,2-diphenylacetaie) and Ci^jjCO^NHS, respectively. TS3 is synthesized by 
reacting T3 with FMOC-Cl (fluorenylmethyloxycarbonylchloride). TS4 is 
synthesized by reacting T3 with tBOCiO (tBOC anhydride or di-t- 

20 butyldicarbonate). TS5, which differs from TSl -4 by having a -H instead of an -I 

at the R'j position, is synthesized by reacting 3,5-diiodothyronine with tBOCjO. 
The general reaction scheme for TSl, TS2, TS3, TS4 and TS5 is depicted in FIG. 
11. It should be noted that in the reaction scheme, both TS5 and its precursor 
both have a hydrogen rather than an iodine at the R', position. 

25 

Synthesis of TS6 and TS7 

TS6 is synthesized by reacting TS5 with paranitrophenylisocyanate. TS7 is 
synthesized by reacting TS6 with TFA (trifluoroaceiic acid), which cleaves the 
tBOC group. These reactions are simple organic synthesis reactions that can be 
30 performed by anyone of ordinary skill m the an. The synthetic scheme for TS6 

and TS7 is diagrammed in FIG. 12. 
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Synthesis of TS8 

TS8 is synthesized by reacting TS5 with Ph.CHNH, (diphenylmethylamine) 
m the presence of tnethylamine and any amide foniimg condensmg reagent, such 
as TBTU (hydroxybenztriazoleuronium tetrafluoroborate) or HBTU 
(hydroxybenztriazoleuromum hexafluorophosphate). The synthesis scheme for 
TS8 is depicted in FIG. 13. 



10 



SYNTHESIS OF 3, 5-011000-3' ISOPROPYLTHVRONINE DERIVATIVES 

For designing a class of antagonists, it is imponant to have a hydrophobic 
group at the 3' position as well as an extension at the 5' position. Preferred 
hydrophobic groups at the 3' position include: methyl, benzyl, phenyl, lodo, and 
heterocyclic stmctures. The synthesis of a 3.5-diiodo-3'-isopropyl-5'-substituted 
thyronine is described below. The example provided describes the specific steps 
for synthesizing the TSIO compound, but this general reaction scheme can be used 
15 by one of ordinary skill in the art to synthesize any number of 3.5,-diiodo-3'- 

isopropyl-5'-substituted thyronine derivatives, which are characterized by havmg 
an extension at the 5' position. Additional compounds of this class can be 
synthesized using known organic synthesis techniques. 

The synthesis of TSIO is described below and is depicted in FIG. 14. 
20 Numbers used in the reaction scheme for TSIO indicating the reaction product for 

each step are in parentheses. 

2-Formvl-6-iso pro pvlanisole (1) : 2-fonnyl-6-isopropylanisole (10.0 g, 61 
mmol), as made by Casiraghi, et al. JCS Perkin I. 1862 (1980) (incorporated by 
reference), is added dropwise to a suspension of sodium hydride (3.7 g, 153 
25 mmoi) in 50 mL THE and 50 mL of DMF in a round bottom flask. The addition 

generates an exothermic reaction and formation of a gray solid. Methyl iodide 
(26.0 g, 183 mmol) is then added dropwise and the reaction mixmre is stirred at 
room temperature for 5 hours. The reaction mixture is quenched with 20 mL of 
water, then poured into 500 mL of water, and is extracted with ether (2 x 300 
mL). The ether layers are combined, washed with water (5 x 1000 mL), dried 
over magnesium sulfate and concentrated in vacuo to provide 10.2 g (94%) of the 
title compound, with the following 'H NMR (CDCIJ propeaies: d 10.30 (s, IH). 
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7.63 (d, IH, J = 3 Hz). 7.50 (d, IH, J = 3 Hz), 7.13 (t. IH. J = 3 Hz). 3.81 (s, 
3H), 3.31 (hepiet. IH. J = 7.5 Hz). 1.19 (d. 6H. J = 7.5 Hz). 

2-f2-Hvdroxvnonvl)-6-isopropvlani$ote (not shown in scheme): 
Octylmagnesium chloride (8.4 mL. 16.9 mmol, 2.0 M) is added dropwise lo a 
solution of 1 (K5 g, 8.4 mmol) in 10 mL THF at -78° C. The reaction mixture is 
stirred for 2 hours with warming to room temperature. The reaction mixture is 
diluted with 50 mL ether and poured into 50 mL water. The ether layer is washed 
with brine (1 x 50 mL), dried over sodium sulfate, and concentrated in vacuo. 
Flash chromatography (silica gel, 10% ether/hexane 15% ether/hexane) 
provides 734 mg (30%) of the title compound with the following 'H NMR (CDCl,) 
properties: d 7.33-7.10 (m, 3H), 5.00 (br. s, IH), 3.81 (s, 3H), 3.33 (heptet, 
IH, J = 7 Hz) 1.90-1-19 (m. 14H), 0.86 (t, 3H, J = 6.5 Hz); HRMS (EI), found: 
292.2404; calc'd: 292.2402. 

2-nonvl-6-isopropvlanisole (2): Compound 2 (663 mg. 2.3 mmol) is 
dissolved in solution of 5 mL ethanol and 5 mL acetic acid, and a spamla tip of 
palladium on carbon catalyst is added. The reaction mixnire is then charged with 
hydrogen gas (using a simple balloon and needle) and the mixture is stirred at 
room temperature overnight. The next day, the reaction mixture is poured into 
ether (100 mL) and the ether layer is extracted with saturated sodium bicarbonate 
(3 x 100 mL) The ether layer is dried over sodium sulfate and concentrated in 
vacuo to provide 581 mg (91%) of (2) with the following 'H NMR (CDCl,) 
properties: d 7.14-7.00 (m, 3H), 3.75 (s, 3H), 3.36 (heptet, IH, J=6.8 Hz), 2.63 
(t, 2H, J = 7,5 Hz), 1.68-1.15 (m. 14H), 0.86 (t, 3H, J = 5.5 Hz); HRMS (EI), 
mass found: 276.2459; calculated: 276.2453. 

Thyronine adduct (4) : Fuming nitric acid (0.071 mL) is added lo 0.184 
mL acetic anhydride chilled to -5'' C. Iodine (66 mg) is added to this mixture 
followed by trifluoroacetic acid (0.124 mL). This mixture is stirred for 1 hour 
with warming to room temperamre, at which point all of the iodine is dissolved. 
The reaction mixmre was then concentrated in vacuo to provide an oily semi-solid 
material. The residue was dissolved in 0.7 mL of acetic anhydride and cooled to - 
20"* C. A solution of anisole (2) (581 mg, 2.1 mmol) in 1.2 mL acetic anhydride 
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and 0.58 mL TFA is added dropwise. The reaction niixmre is stirred at -20= for 
1 hour, then stirred overnight with warming to room temperature. The reaction 
mixture is partitioned between water and methylene chloride. The methylene 
chloride layer is dried over sodium sulfate and concentrated in vacuo to provide 
the iodonium salt (3) as an oil. This material is not purified or characterized, and 
is directly introduced into the coupling reaction. 

N-Trifluoroacetyl-3,5-diiodotyrosine methyl ester (552 mg, 1.0 mmol) 
prepared according to the procedure of N. Uwis and P. Wallbank, Synthesis 1103 
(1987) (incorporated by reference) and all of the crude iodonium salt (3) from 
above is dissolved in 5 mL of anhydrous methanol. Diazabicyclo[5.4.01undecane 
(DBU) (183 mg. L2 mmol) and a spatula tip of copper-bronze are added and the 
resulting mixture is stirred at room temperamre overnight. The next day, the 
reaction mixture is filtered, and the filtrate is concentrated in vacuo. The cmde 
residue is purified by flash chromatography (silica gel, 10% ethyl acetate/hexane) 
to provide 30 mg (4%) of the protected thyronine adduct (4). 

De protected thyronine fTSlO) : The protected thyronine 4 (30 mg. 0.04 
mmol) is dissolved in a mixture of 2.25 mL acetic acid and 2.25 mL 49% 
hydrobromic acid. The reaction mixture is heated to reflux for 5 hours. The 
reaction mixture is cooled to room temperamre, and the solvents are removed in 
vacuo. Water is added to triturate the oily residue into a gray solid. This solid 
material is filtered, washed with water, and dried over P,0, in vacuo to provide 24 
mg (81%) of the title compound, TSIO, with the following 'H NMR (CDCl,) 
propenies: d 7.57 (s, IH), 6.86 (s, IH). 6.45 (s, IH). 6.34 (s, IH). 4.81 (m, 
IH), 3.86 (s, 3H), 3.71 (s, 3H). 3.33-3.05 (m. 3H). 2.58-2.47 (m, 2H). L62-0.76 
(m, 23H); MS (LSIMS): = 817.0. 

As mentioned above, this reaction scheme can be modified by one of 
ordinary skill in the art to synthesize a class of compounds characterized by 3.5- 
diiodo-3'isopropylthyronine derivatives, wherein (1) the 3' isopropyl group can be 
replaced with a hydrophobic group, including methyl, benzyl, phenyl, lodo. and 
heterocyclic structures, and (2) a wide variety of chemical stoictures can be 
incorporated at the 5' position, including alkyl groups, planar aryl, heterocyclic 
groups, or polar and/or charged groups. 
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The aldehyde (1) in the above reaction scheme is a versatile synthetic 
intermediate which allows for the attachment of a variety of chemical moieties to 
the 5' position of the final thyronine derivative. In addition, a variety of chemical 
reactions can be used to attach the chemical moieties. These reactions are well 
known in the art and include organometallic additions to the aldehyde (including 
Grignard reagents, organolithiums, etc.), reductive amination reactions of the 
aldehyde with a primary or secondary amine, and Wittig olefination reactions with 
a phosphorous ylid or stabilized phosphonate anion. Other possibilities include 
reduction of the aldehyde to a benzyl alcohol allowing for etherification reactions 
at the 5' position. As mentioned above, these methods allow for a wide variety of 
chemical structures to be incorporated at the 5' position of the fmal thyronine 
derivative, including alkyl groups, planar aryl, heterocyclic groups or polar and/or 
charged groups. 

Synthesis of 3, 5-dibromo-4-(3',5'-diisopropyl-4'-hydroxyphenoxy) benzoic acid 
(Compound 11). 



(a) A mixture of 2,6-diisopropyl phenol (20 g, 0.11 mol), potassium 
carbonate (62 g, 0.45 mol), acetone (160 ml) and methyl iodide (28 ml, 0.45 
mole) is refluxed for three days. The reaction mixmre is filtered through celite, 
evaporated, dissolved in ether, washed twice with IM sodium hydroxide, dried 



2,6-diisopropyl anisole as a slightly yellow oil. 

(b) Fuming nitric acid (12.4 ml. 265 mmol) is added dropwise to 31.4 ml 
of acetic anhydride which is cooled in a dry ice/carbon tetrachloride bath. Iodine 




over magnesium sulphate and concentrated to afford 15.1 g (0.08 mol, 70%) of 



97/21993 PCT/US96/20778 

49. 

11.3 g. 44.4 mmol) is added in one ponion followed by dropwise addition of 
trifluoroacetic acid (20.5 ml, 266 nimole). The reaction mixture is stirred at room 
temperature until all the iodme is dissolved. Nitrogen oxides are removed by 
flushing nitrogen into the vessel. The reaction mixture is concentrated, the residue 
is dissolved in 126 ml of acetic anhydride and is cooled m a dry ice/carbon 
tetrachloride bath. To the stirred solution 2,6-diisopropylanisole (51 g, 266 mmol) 
in 150 ml of acetic anhydride and 22.6 ml of trifluoroacetic acid is added 
dropwise. The reaction mixture is left to stand at room temperature over night 
and then is concentrated. The residue is taken up in 150 ml of methanol and 
treated with 150 ml of 10% aqueous sodium bisulfite solution and 1 liter of 2M 
sodium borotetrafluoride solution. After the precipitate aggregates, petroleum 
ether is added and the supernatant is decanted. The precipitate is triturated with 
petroleum ether, filtered, washed with petroleum ether and dried at room 
temperature in vacuo. This affords 34 g (57 mmol, 65%) of bis(3,5-diisopropyl-4- 
methoxyphenyDiodonium tetrafluoroborate as a white solid. 

(c) To a stirred solution of 3,5-dibromo-4-hydroxybenzoic acid (12 g, 40,5 
mmol) in 250 ml of methanol, thionyl chloride (3 ml) is added dropwise. The 
reaction mixmre is refluxed for five days, water is added and the precipitated 
product is filtered off. The residue is dissolved in ethyl acetate. From the 
aqueous phase, methanol is removed by concentration. The aqueous phase is then 
samrated with sodium chloride, and extracted with ethyl acetate. The combined 
organic phases are dried over magnesium sulphate, filtered and concentrated. This 
gives 12.5 g (40.5 mmol, 100%) of 3,5-dibromo-4-hydroxymethyl benzoate as a 

white crystalline solid. 

(d) The products obtained in steps b and c are reacted with each other 
according to the following protocol. To bis(3.5-diisopropyM- 
methoxyphenyl)iodonium tetrafluoroborate (2.86 g. 4.8 mmoie) and copper bronze 
(0.42 g. 6.4 mmole) in 7 ml. of dichloromethane at 0° C is added dropwise a 
solution of 3,5-dibromo-4-hydroxymethyl benzoate (1.0 g. 3.2 mmole) and 
triethylamine (0.36 g. 3.5 mmole) in 5 ml of dichloromethane. The reaction 
mixmre is stirred in the dark for eight days and then is filtered through celite. 
The nitrate is concentrated and the residue is purified by column chromatography 
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(silica gel, 97:3 petroleum ether/ethyl acetate) to give 0.62 g (1.2 mmole. 39%) of 
3,5-dibromo-4-(3'.5'-ciiisopropyl-4'-methoxyphenoxy)methyl benzoate as a solid. 

(e) The product from step d (0.2 g, 0.4 mmole) Ls dissolved in 2 ml. 
dichloromethane, is put under nitrogen and is cooled at -40° C. To the stirred 
solution is added IM BBr, (L2 ml. 1.2 mmole) dropwise. The reaction mixture is 
allowed to reach room temperature and then is left over night. It is cooled to 0° 
C and then hydrolyzed with water. Dichloromethane is removed by concentration 
and the aqueous phase is extracted with ethyl acetate. The organic phase is 
washed with IM hydrochloric acid and brine. Then it is dried over magnesium 
sulphate, filtered and concentrated. The residue is chromatographed (silica, 
96:3.6:0.4 dichloromethane/methanol/acetic acid) producing 93 mg (0.2 mmole, 
51%) of 3,5-dibromo-4-(3',5'-diisopropyI-4'-hydroxyphenoxy)benzoic acid as a 



white solid. nmr (CDClj) 5 1.23 (d, 12H, methyl), 3.11 (m, 2H, CH), 6.50 
(s. 2H, 2.6-H) 8.33 (s, 2H. 2',6'-H). 

TABLE 1 and FIG. 15 depict the structures of several TR ligands. 
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Example 2 - Receptor Binding Assays of TR Ligands 

To test the ability of synthesized TR ligands to bind to a thyroid receptor (TR), 
the binding affinity of a TR ligand for TR is assayed using TR's prepared from rat liver 
nuclei and 125, T, as described in J.D. Apriletti, J.B. Baxter, and T.N. Lavin, 7. Biol. 
5 Chem., 263: 9409-9417 (1988). The apparent Kd's are calculated using the method 
described by Apriletti (1995) and Apriletti (1988). The apparent Kd's are presented in 
TABLE 2. The apparent Kd's (App.Kd) are determined in the presence of the sample to 
be assayed. 1 nM ['"IJTj, and 50/ig/ml core histones, in buffer E (4(X) mM KCI, 200 
mM potassium phosphate, pH 8.0, 0.5 mM EDTA, 1 mM MgCl^, 10% glycerol, 1 mM 
10 DTT) in a volume of 0.21 ml. After incubation overnight at 4°C, 0.2 ml of the 

incubation mixmre is loaded onto a Quick-Sep Sephadex G-25 column (2,7 x 0.9 cm, 1.7 
ml bed volume) equilibrated with buffer E. The excluded peak of protein-bound ['^IJTj 
is eluted with 1 ml of buffer E, collected in a test tube, and counted. Specific T3 bindmg 
is calculated by subtracting nonspecific binding from total binding. 



15 
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20 Example 3 - Increased Nuclear Protein Coactivation by TR Ligands 

To test the ability of TR ligands to activate the binding of TR to the nuclear 

activation protein RIP-140 (a nuclear protein that can bind to nuclear receptors, such a; 

the estrogen receptor), a TR ligand is liganded to TR and then incubated wtth RIP-140 

described in V. Cavailles, et ai. EMBO J.. 14(15):3741- 3751 (1995), which is 
25 incorporated by reference herein. In this assay, 35s-RIP-140 protein binds to liganded 

but not unliganded TR. Many TR 35s ligands can activate RIP-140 binding as shown i 

TABLE 2. 
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Example 4 - TR ligand binding and TR activation i> Cultl-red Cells 

To test TR activation of transcription in a cellular environment, TR ligands are 
assayed for their ability to activate a reporter gene, chloramphenicol transferase ("CAT"), 
which has a TR DNA binding sequence operatively linked to it. Either GC or L937 cells 

5 (available from the ATCC) can be used, respectively). In such assays, a TR ligand 
crosses the cell membrane, binds to the TR, and activates the TR, which in turn activates 
gene transcription of the CAT by binding the TR DNA binding region upstream of the 
CAT gene. The effective concentration for half maximal gene activation (EC50) is 
determined by assaying CAT gene activation at various concentrations as described herem 

10 and in the literamre. The results of CAT gene activation experiments are shown in 
TABLE 2 

CAT GENE ACTIVATION ASSAYS 

Functional response to thyroid hormone (3»5,3'-triiodo-L-thyronine, T3) and TR 
15 ligands is assessed either in a rat pituitary cell line, GC cells, that contain endogenous 
thyroid hormone receptors (TRs) or U937 cells that contain exogenous TRs expressed as 
known in the an. GC cells are grown in 10-cm dishes in RPMI 1640 with 10% newborn 
bovine serum, 2 mM gluiamine, 50 units/ml penicillin and 50 ^g/m\ streptomycin. For 
transfeciions, cells are trypsinized, resuspended in buffer (PBS, 0.1% glucose) and mixed 
20 with a TREtkCAT plasmid (10 mg) or phage in 0.5 ml buffer (15 ±5 million cells) and 
electroporated using a Bio-Rad gene pulser at 0.33 kvolts and 960 mF. The TREtkCAT 
plasmid contains two copies of a T, response element (AGGTCAcaggAGGTCA) cloned in 
the Hind III site of the pUC19 polylinker immediately upstream of a minimal (-32/+45) 
thymidine kinase promoter linked to CAT (tkCAT) coding sequences. After 
25 electroporation, cells are pooled in growth medium (RPMI with 10% charcoal -treated, 
hormone stripped, newborn bovine serum), plated in 6-well dishes and treated with either 
ethanol or hormone. CAT activity is determined 24 hours later as described D. C. 
Leitman, R. C. J. Ribeiro. E. R. Mackow, J. D. Baxter, B. L. West, J. Biol. Chem. 
266, 9343 (1991), which is incorporated by reference herein. 



30 
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Effect of TS-10 on the Transcriptional Reglxation of the DR4-ALP 
Reporter Gene in the Presence or Absence of T3. 

Characteristics of the TRAF ceils: TRAFal are CHO Kl cells stably transformed with 
an expression vector encoding the human thyroid hormone receptor a 1 and a DR4.ALP 
5 reporter vector; TRAFbl are CHO Kl cells stably transformed with an expression vector 
encodmg the human thyroid hormone receptor jSl and a DR4-ALP reporter vector. 

Interpretation of the effect of compound TS-10 on the transcriptional regulation of 
the DR4-ALP reporter gene in the presence or absence of T3. 

10 

TRAFal reporter cells: TS-10 alone (open circles) induces a partial activation of the 
expression of the ALP reporter protein amounting to approximately 27 % of the maximal 
effect by the natural thyroid hormone T3. In the presence of T3 (filled circles). TS-10 
has a weak antagonistic effect. The EC50 concentration for the agonistic effect of TS-10 
15 and the EC50 concentration for its T3 antagonistic effect, respectively, is indicated in 
FIG. 18 

In FIG. 18, open and filled circles with dotted lines show the dose-dependent effect of 
TS-10/T3 on the toxicity marker (MTS/PMS). reduction of tetrazolium salt in the 
20 mitochondria, displayed on the right y-axis as optical density. There is no obvious toxic 
effect of TS-10 on the MTS-PMS marker but there is a clear effect on the morphology of 
the cells, as can be seen under the light microscope, at the highest concentration of TS-10 
(32 mM) both in the absence and presence of T3, respectively (not shown in the figure). 

25 TRAFbl reporter cells: TS-10 alone (open circles) induces a partial activation of the 
expression of the ALP reporter protein amounting to approximately 35% of the maximal 
effect by T3. The EC50 concentration for the agonistic effect of TS-10 is indicated in 
FIG. 19. In the presence of T3 (filled circles). TS-10 shows, if anything, a slight 
potentiation of the T3 effect on the expression of the ALP reporter protein. The T3 

30 inhibitory effect of TS-10 at its highest concentration used (32 mM) is a toxic effect 
rather than T3 antagonism. 
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In FIG. 19, open and filled circles with dotted lines show the dose-dependent effect of 
TS-10/T3 on the toxicity marker (MTS/PMS), reduction of tetrazolium salt in the 
mitochondria, displayed on the right y-axis as optical density. There is no obvious toxic 
effect of TS-10 on the MTS-PMS marker but a clear effect on the morphology of the cells 
5 can be observed, under the light microscope, at the highest concentration of TS-10 (32 
mM) both in the absence and presence of T3, respectively (not shown in the figure). 

HepG2 (HAF18) reporter cells: TS-10 alone (open circles) induces a partial activation 
of the expression of the ALP reponer protein amounting to slightly more than 50% of the 

10 maximal effect by T3. The EC50 concentration for the agonistic effect of TS-10 is 
indicated in FIG, 20. In the presence of T3 (filled circles), TS-10 shows no effect i.e. 
no T3 antagonism nor potentiation/additive effect to T3. Open and filled circles with 
dotted lines show the dose-dependent effect of TS-10/T3 on the toxicity marker 
(MTS/PMS), reduction of tetrazolium salt in the mitochondria, displayed on the right y- 

15 axis as optical density. There is no obvious toxic effect of TS-10 on the MTS/PMS 

marker or on the morphology of the cells, as can be observed using a light microscope, at 
any concentration of TS-10/T3 used. 

Example 5 - Comparisons of Human TR-of and Human TR-^ 
20 Competition for [^"qTj binding to TR LBD by and Triac 

The drug, triac. is a thyroid hormone agonist. Triac is 3,5.3'-triiodothyroacetic 

acid and is described in Jorgensen, Thyroid Hormones and Analogs in 6 Hormonal 

Proteins and Peptides. Thyroid Hormones at 150-151 (1978). Another compound that can 

be used in place of triac is 3,5-diiodo-3'-isopropylthyroacetic acid. Competition assays 
25 are performed to compare the displacement of [^^I]T, from binding with human TR-a 

LBD or human TR-i3 LBD by unlabeled T, or triac. The results of such assays are 

depicted in FIG. 16. 

Standard binding reactions are prepared contaimng 1 nM [^"IjT^, 30 fmol of 

human TR-a (empty symbols) or /? (solid symbols), and various concentrations of 
30 competing unlabeled T, (circles) or triac (triangles). Assays are performed m duplicate. 
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Scatchard Analysis of V''l]T, Binding to TR 

Human TR-a (left panel) or human TR-^ (right panel) is assayed for Tj bindmg m 
the presence of increasing concentrations of ['^^1X3. The apparent equilibrium 
dissociation constant (20 pM for and 67 pM for 0) is calculated by Imear regression 

5 analysis and is depicted in FIG. 17. 

3, 5-DlBROMO-4-(3\5^D[ISOPROPYL-4'-HYDROXYPHENOXY) BENZOIC ACID IS A 

TR-A Selective SvimiETic Ligand. 



10 




HO 



15 3, 5-dibromo-4-(3\5'-diisopropy!-4'-hydroxyphenoxy) benzoic acid (Compound 

11), the strucnire of which is drawn above, is assayed for bindmg to the two different 
isoforms of the TR. Tra and TR^. Compound U exhibits an IC50 of 1.6 /xM for 
binding to TRa and an IC50 of 0.91 ^iM for binding to TR/?. Assays for determining 
selective binding to the TRa or TRj3 LBD can include reporter assays, as described 

20 herem. See also Hollenberg, et aL, J. Biol. Chem., 270(24)14274-14280 (1995). 



Example 6 - Preparation and Purification of a TR-a LBD 

Rat TR-a LBD, residues Metl22 - Val4l0, is purified from £. coli ("LBD- 
122/410"). The expression vector encoding the rat TR-a LBD is freshly transfected mto 

25 E. coli strain BL21(DE3) and grown at 22X in a 50-titer fermemer using 2x LB 

medium. At an of 2.5-3, IPTG is added to 0.5 mM and growth is continued for 3 h 
before harvesting. The bacterial pellet is quickly frozen in liquid nitrogen and stored at - 
70"C until processed. Extraction and purification steps are carried out at 4°C. The 
bacteria are thawed in extraction buffer (20MM Hepes, pH 8.-. 1 mM EDTA, 0.1% 

30 MTG. 0.1 mM PMSF. and 10% glycerol) at a ratio of 10 ml buffer/g bacteria. Bacteria 
are lysed by incubation for 15 min. with 0.2 mg/ml lysozyme and sonicated at maximum 
power while simultaneously homogenized with a Brinkmann homogenizer (Model PT 
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10/35 with generator PTA 35/2) until the solution loses its viscosity. After centrifugation 
for 10 min at 10,000 g. the supernatant is adjusted to 0.4 M KC!, treated with 0.6% PEI 
to precipitate fragmented DNA, and centrifuged for 10 min at 10,000 g. The rat TR-a 
LBD in the supernatant is then precipitated with 50% ammonium sulfate and centrifuged 
5 for 10 min at 10,000 g. The precipitate is resuspended with buffer B (20 mAl Hepes. pH 
8.0, 1 mM EDTA, 1 mM DTT. 0.1 mM PMSF, 0.01% Lubrol, and 10% glycerol) to a 
final conductivity of 9 mS/cm (approx. 0.7 M ammonium sulfate) and centrifuged 1 h at 
100,000g. The supernatant is frozen in liquid nitrogen and stored at -70°C. 

The crude extract is thawed, bound with a tracer amount of [^'"I]T-,, and loaded 

10 directly onto a phenyl-Toyopearl hydrophobic interaction column (2.6 x 18 cm, 95 ml bed 
volume) at 1.5 ml/min. The column is eluted with a 2-h gradient from 0.7 ammonium 
sulfate, no glycerol to no salt, 20% glycerol in buffer C (20 mM Hepes, pH 8.0, 0.5 mM 
EDTA, 1 mM DTT, 0.2 mM PMSF). The rat TR-a LBD prebound to tracer ['"IJT, 
(less than 0.005% of total rat TR-a LBD) is detected using a flow-through gamma 

15 emission detector, whereas unliganded rat TR-a LBD is assayed by postcolumn ['"rjTj 
binding assays (described herein). 

The phenyl-Toyopearl unliganded rat TR-a LBD peak fractions are pooled, diluted 
with buffer B to a conductivity of 0.5 mS/cm (equivalent to approx. 20 mM ammonium 
sulfate), loaded onto a TSK-DEAE anion-exchange column (2x15 cm, 47 ml bed 

20 volume) at 4 ml/min, and eluted with a 60-min gradient from 50 to 200 mM NaCI in 
buffer B. 

The unliganded rat TR-a LBD peak fractions from TSK-DEAE are pooled, diluted 
twofold with buffer B, loaded at 0.75 ml/min on a TSK-heparin HPLC column (0.8 x 7.5 
cm, 3 ml bed volume), and eluted with a 50 to 400 mM NaCl gradient in buffer B. 

25 The pool of unliganded rat TR-a LBD peak fractions from the TSK-heparin 

column is adjusted to 0.7 M ammonium sulfate, loaded at 0.75 ml/min on a TSK-phenyl 
HPLC column (0.8 x 7.5 cm, 3 ml bed volume), and eluted with a 60-min gradient from 
0.7 M ammonium sulfate without glycerol to no salt with 20% glycerol in buffer C The 
fractions containing unliganded rat TR-a LBD are pooled and incubated with a five fold 

30 excess of hormone for 1 h, the salt concentration is adjusted to 0.7 M ammonium sulfate, 
and the sample is reloaded and chromatographed on the same column as described above. 
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EXAMPLE 7 - Crystallization of Liganded TR-a LBD 

Material from a single LBD-122/410 preparation is divided into batches, and 
quantitatively bound with one of the following ligands: Dimii. T,. or triac IpBr, 
(3,5dibromo-3'isopropylthyronine) for the final purification step. 
5 To maintain full saturation of rat TR-a LBD with a iigand, and to prepare the 

complex for crystallization, the ligand-bound rat TR-a LBD is concentrated and desalted 
in an Amicon Centricon-10 microconcentrator (McGrath et al, Biotechniques , 7:246-247 
(1989), incorporated by reference herein), using 10 mM Hepes (pH 7.0), 3.0 naM DTT. 
and 1.0 nM to 10 nM ligand. 
10 Factorial crystallization screening trials (Jancarik & Kim, J. Appl. Crystallogr. 

24:409-411 (1991) incorporated by reference herein) are carried out for rat TR-a LBD 
bound to selected ligands using hanging-drop vapor diffusion at 1 7^*0 (with 1 fi\ protein 
solution, 1 ti\ precipitant solution and a 0.5 ml reservoir using silanized coverslip: 
(McPherson, Preparation and Analysis of Protein Crystals (1982), incorporated by 
15 reference herein). Rat TR-a LBD is not stable at 4°C and is stored at -80**C, where it 
maintains its avidity for hormone and its crystallizability for approximately two to three 
months. These procedures are carried out as described in McGrath, M.E. et al., J. MoL 
Biol. 237:236-239 (1994) (incorporated by reference).) Crystals are obtained in condition 
21 of the screening trials (Jancarik & Kim 1991) and conditions are then optimized. 
20 Wedge-shaped crystals are reproducibly obtained with hanging-drop vapor fusion at 22°C 
with 15% 2-methyl-2,4-pentanediol (MPD), 0.2 M ammonium acetate and 0.1 M sodium 
cacodylate (pH 6.7), 3 mM DTT, with 2 ^\ protein solution, 1 ^1 precipitant solution and 
a 0.6 ml reservoir using silanized coverslip, and with 8.7 mg/ml (Dimit), 5.5 mg/ml 
(IpBr^), 5 mg/ml (triac), or 2.3 mg/ml (T3) over a period of three days. Under these 
25 conditions, diffraction quality crystals (dmiension 0.5 x 0.2 x 0.0075 mm^) can be grown 
at ambient temperature (22°C). The best crystals have a limiting dimension of 
approximately 100 ^m and are obtained at a protein concentration between 2.3 and 8.7 
mg/ml in the presence of 3 mM DTT. The crystals are of the monoclinic space group 
C2, with one monomer in the asymmetric unit. 

30 
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Example 8 - Crystallization of Human TR-jS LBD Complexed with T3 or 
Triac 

Human TR-/3 LBD complexed with T, and human TR-p LBD complexed with 
triac are purified according to the same procedures described above for the rat TR-a 
5 LBD, with the following modifications. 

The expression of human TR-jS LBD differs from the rat TR-a LBD in that the 
human TR-^ LBD residues extend from the amino acid at position 716 through the amino 
acid at position 1022, according to the amino acid numbering scheme for the various 
nuclear receptor LBDs depicted in FIG. 3. FIG. 3 illustrates a numbering scheme 

10 applicable to all of the nuclear receptors listed as well as to any additional homologous 
nuclear receptors. The vertical lines on FIG. 3 at position 725 and at position 1025 
delineate the preferred minimum amino acid sequence necessary to obtain adequate 
binding of ligand. The amino acid sequence from position 716 to position 1022 according 
to the numbering scheme of FIG. 3 corresponds to the amino acid positions 202 to 461 

15 according to the conventional numbering of the amino acid sequence of human TR-/3 
which is publicly available. Also, the human TR-jS LBD is expressed with a histidine 
tag, as described in Crowe et al.. Methods in Molecular Biology 31:371-387 (1994), 
incorporated by reference herein. 

The purification of human TR-^ LBD is the same as that described above for the 

20 rat TR-a LBD with the following exceptions. First, before the purification step using the 
hydrophobic interaction column, a step is added in which the expressed human TR-^ LBD 
is purified using a nickel NTA column (commercially available from Qiagen, Chatsworth, 
CA) according to manufacturer's instructions, and eluted with 200 mM imidazole. The 
second difference is that in the purification of the human TR-/? LBD, the purification step 

25 using a heparin column is omitted. 

The crystallization of human TR-/3 LBD bound to T3 or triac is as follows. 
Crystals are obtained in condition 7 of the factorial screen using hanging drops as before 
at ambient temperature (22'C) using the factorial crystallization screening trials of 
Jancarik &. Kim (1991) and using the commercially available product from Hampton 

30 Research, Riverside). The following are optimum conditions: hexagonal bipyrimidal 
crystals are grown at 4''C for 2-3 days from hanging drops containing 1.0-1.2 M sodium 
acetate (pH unadjusted) and 0.1 M sodium cacodylate (pH 7.4), 3 mM DTT, with either 
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a I protein solution, 1 m1 precipitant solution or 2 ^xl protein solution, 1 ^1 precipitant 
solution and a 0.6 ml reservoir using silanized coverslip, at a protein concentration of 1- 
10 mg/ml. The best crystals have a limiting dimension of 200 ^im. 

The crystal system for human TR-/? LBD bound to Tj or triac is trigonal with the 
5 space group p3,2i. The umt cell dimensions are cell length a = cell length b = 68.448 
angstroms, cell length c = 130.559 angstroms. The angles are a = 90\ (3 = 90°, 
gamma - 120°. 

Example 9 - Determination of Liganded TR-a LBD Crystal Structure 

10 Data from each of three cocrystals (Rat TR-a LBD with Dimit, T3 and IpBrj) is 

measured on a Mar area detector at Stanford Synchrotron Radiation Laboratory beamline 
7_1 (X = 1.08 angstroms) using 1.2' oscillations. 

Data from the Tj cocrystal is measured with the b* axis approximately parallel 
with the spindle. The crystals are flash frozen at -178°C in a nitrogen gas stream with the 

15 MPD mother liquor serving as the cryosolveni. An orientation matrix for each crystal is 
determined using REFIX (Kabsch. W., J. Appi Crysiallogr. 26:795-800 (1993) 
incorporated by reference). Reflections are integrated with DENZO (commercially 
available from Molecular Structure Corp., The Woodlands, Texas), and are scaled with 
SCALEPACK (as described in Otwinowski, Z. Proceedings of (he CCP4 Study Weekend: 

20 "Data Collection and Processing, " 56-62 (SERC Daresbury Laboratory, Warrington, UK 
1993) incorporated by reference). 

For the T^ data set, Bijvoet pairs are kept separate, and are locally scaled using 
MADSYS (W. Hendrickson (Columbia University) and W. Weis (Stanford University)). 
Cocrystals prepared from the three isosteric ligands are isomorphous. MIR 

25 analysis is performed using programs from the CCP4 suite (Collaborative Computational 
Project, N.R. Acta Crystallogr. D50:760-763 (1994), incorporated by reference herein). 
Difference Pattersons is calculated for both Tj and IpBr., taking the Dimit cocrystal as the 
parent. The positions of the three iodine atoms in the T, difference Patterson are 
unambiguously determined from the Marker section of the density map as peaks of 11a 

30 above background. The positions for the two bromine atoms in the IpBr, cocrystals, are 
located independently, as peaks 8a above the noise level. Phases for the LBD-122y'410 are 
calculated from the solution to the IpBr. difference Patterson, and are used to confirm the 
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location of the unique third iodine of the T, cocrystal. Halogen positions are refined with 
MLPHARE. including the anomalous contributions from the iodine atoms (Otwmowski, 
Z, Proceedings of the CCPR Study Weekend 80-86 (SERC Daresbury Laboratoo'. 
Warrington, UK 1991)). The MIRAS phases are improved through solvent 
5 flattening/histogram matching using DM (Cowtan, K., Joint CCP4 and ESF EACBM 
Newsletter on Protein Crystallography 31: 34-38 (1994), incorporated by reference 
herein). 

A model of the LBD- 122/4 10 with Dimit bound is built with the program O from 
the solvent flattened MIRAS 2.5 angstrom electron density map (Jones et al.. Acta 

10 Crystallogr. A 47:110-119 (1991), incorporated by reference herein). The initial model, 
without ligand, (Rcryst = 40.1%), is refined using least-squares protocols with XPLOR. 
The Dimit ligand is built into unambiguous Fo-Fc difference density during the following 
round. Subsequent refinement employs both least-squares and simulated annealing 
protocols with XPLOR (Brunger et al.. Science 235:458-460 (1987), incorporated by 

15 reference herein). Individual atomic B-factors are refined isotropically. As defined in 
PROCHECK, all residues are in allowed main-chain torsion angle regions as described in 
Laskowski et al., J. Appl. Crystallogr. 26:283-291 (1993). incorporated by reference 
herein. The current model is missing 34 residues (Meti^^-Glni,^) at the N-tenninus, and 5 
residues (Glu406-VaIj,o) at the C-terminus. 

20 In addition, the following residues are not modeled beyond due to poor 



density': 184, 186, 190. 198, 206, 209, 240, 301, 330, 337, 340. 343, 359, and 395. 



122/410, residues Argij^-SeriBj. Trpigj-Glyig^. Serj^-Asp^o^ and Asp3o8-Phe405; three 
cacodylate-modified cysteines: Cysjjj. Cysjgo and Cysj,.; and 73 solvent molecules 
25 modeled as water (2003 atoms). 

*R,,, = 100 X EsM I. I I.- 11/ E^u E.I. 
tR^, - 100 X L,JFp„-FkI/ 

The occupancy for the two bromine sites is set to 35 electrons. The occupancies of the 
iodine sites are relative to this value. 
30 §Phasing power = (FH) / (t), where (FH> is the mean calculated heavy atom structure 
factor amplitude and (e) is the mean estimated lack of closure. 



The average B- value for protein atoms is 34.5 A'. The final model consists of the LBD- 
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llRcuIlis = (e> / (iso), where (e) is the mean estimated lack of closure and (iso) is the 
isomorphous difference. 

IRcryst = 100 x E,,, IF,-Fd / IFJ where F„ and R are the observed and 
calculated structure factor amplimdes (for data F/a > 2). The Rfree was calculated usmg 
5 3% of the data, chosen randomly, and omitted from the refmement. 

§ Correlation coefficient = (IFj - !F,[) x (IFj - [Fj)/i:,„ (tFj - lFj)^xI,„ 

EXAMPLE 10. Phasing of the rTRa LBD complex with Triac 

10 Due to the possible non-isomorphism of the rTRa LBD complex with Triac, a molecular 
replacement solution is determined using AMORE (Navaza, J., Acta Crystal log raphica 
Section A-Fundamentals of Crystallography 50:157-63 (1994) from a starting model 
consisting of rTRa LBD complex with T^, but with the ligand, all water molecules, and 
the following residues omitted: Asn 179. Arg228, Arg262, Arg266, and Ser 277. Strong 

15 peaks are obtained in both the rotation and translation searches, with no significant (> 
0,5 times the top peak) false solutions observed (Table 3). Strong positive density present 
in both the anomalous and conventional difference Fourier maps confirm the solution. 
Maps are calculated using sigma-A weighted coefficients output by REFMAC 
(Murshudov. et al. "Application of Maximum Likelihood Refinements," in the 

20 Refinement of Protein Strucmres, Proceedings of Daresbury Saidy Weekend (1996)) after 
15 cycles of maximum likelihood refmement. Triac, the omitted residues, and water 
molecules 503, 504, 534 (following the numbering convention for the TR complex with 
T3) are built into the resulting difference density using O (Jones et. al ); the 
confonnations of these residues are further confirmed in a simulated-annealing omit map 

25 (Brunger et. al.). The complete model is then refined using positional least-squares, 
simulated annealing, and restrained, grouped B factor refinement in XPLOR to an Rcryst 
of 23.6% and an Rfree of 24.1% 

Example 11. Connecting QSAR with Structure in the Thyroid HoRiMONE 
30 Receptor 

The conclusions of classic thyroid hormone receptor quantitative strucmre-aciivity 
relationships may be summarized as follows: 
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1) the R4'-hyciroxyl group functions as a hydrogen bond donor, 

2) the amino-propionic acid interacts electrostatically through the carboxylate 
anion with a positively charged residue from the receptor; 

3) the preferences of Rj/Rs substiment are I>Br>Me> >H; 

5 4) the preferences of the Rj'-substituem are Ipr> I>Br>Me> >H. 

The structure of the thyroid hormone receptor ligand binding domam complexed with the 
agonists 3,5,3*-triiodothyronine (T3). 3,5-dibromo-3'-isopropylthyromne (IpBr.), 3,5- 
dimethyI-3'MSopropylthyronine (Dimit), and 3,5.3'-triiodothyroacetic acid (Triac), as 

provided herein, permits: 
10 1) the identification of receptor determinants of binding at the level of the 

hydrogen bond; 

2) the association of these determinants with the predictions of classic thyroid 
hormone receptor QSAR; and 

3) prediction as to which determinants of binding are rigid, and which are 
15 flexible, for both the ligand and the receptor. 

This classification for the agonists of the type (Ri = amino-propionic, acetic acid; 
R3,Rs = I,Br,Me; Rj'^Ipr,!) is given below (for the representative ligand T,); 



F = Fiducial (always satisfied) 
20 A = Adjustable 




30 
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Based upon the methods and data described herein, the following is an 
embodiment of the computational methods of the invention, which permit design of 
nuclear receptor ligands based upon interactions between the structure of the amino acid 
residues of the receptor LBD and the four different ligands described herein. The small 
5 molecule stmctures for the ligands can be obtained from Cambridge Structural Database 
(CSD), and three dimensional models can be constructed using the methods described 
throughout the specification. The following are factors to consider in designing synthetic 
ligands: 

1) Histidine 381 acts as a hydrogen bond acceptor for the R,' hydroxy 1, with the 

10 optimal uutomer maintained by water molecules. See FIG. 23 and FIG 24. Histidine is 
the only hydrophilic residue in this hydrophobic pocket that surrounds the R,' substituent. 
Histidine can be either a hydrogen bond acceptor or donor, depending on its tautomeric 
state. It is preferably a hydrogen bond donor, but can tolerate being a hydrogen bond 
acceptor, as for example, when there is a methoxy at the R,' position of the ligand; 

15 2) Arginines 228, 262, and 266 interact directly and through water-mediated 

hydrogen bonds with the R, -substituent, with the electrostatic interaction provided by 
Arginine 266 (as in the Triac complex). This polar pocket is illustrated by FIG. 23 - 
FIG. 25. FIG. 23 depicts in the TRa ligand binding cavity, where T3's amino- 
propionic Rl- substituent interacts with Arg 228, H0H502. H9H503 and H0H5(>4 via 

20 hydrogen bonds. FIG. 24 depicts triac in the ligand binding cavity, with its -COOH R, 
substituent in the polar pocket. In FIG. 24. Arg 228 no longer shares a hydrogen bond 
with the ligand, but the -COOH substiaient forms hydrogen bonds with Arg 266. 
FIG. 25 superimposes and triac in the ligand binding cavity and shows several 
positionally unchanged amino acids and water molecules, and selected changed interacting 

25 amino acids and water molecules. The three figures illustrate pans of the polar pocket 
that can change and those parts that do not move upon binding of different ligands. For 
example, the Arg 262 at the top of the polar pocket does not move, even when the R, 
substituent has changed from a -COOH to an aminopropionic acid group. However, the 
other two Arginines. Arg 228 and Arg 266. demonstrate nexibility in the polar pocket to 

30 respond to the change in the size or chemical naure of the R, subsiiment.; 

3) Inner and outer pockets for the R3/R5 substituems are fonned by Ser260, 
Ala263, lle299; and Phe 218, lle221, lle222. respectively. See FIGS. 21 and 22. The 
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inner pocket is filled by either the R3 or the R5 substiment. regardless of the size of the 
substituent, and may act as a binding determinant by positioning the hgand in the 
receptor. Optimally, the inner pocket amino acids interact with an R3 or R5 substiment 
that is no larger than an iodo group. If the inner pocket is filled by the R3 substituent, 

5 then the outer pocket interacts with the R5 substiment and vice versa. The outer pocket 
can adjust to the size of its substiment through main chain motion centered at the break in 
helix 3 (Lys220-lle221). suggesting that the bending of H3, and motion of the N-terminal 
portion of H3, may represent a conformational change induced on ligand binding. The 
outer pocket has greater flexibility than does the inner pocket in terms of accommodating 

10 a larger substiment group. 

4) A pocket for the R3'-substiment is formed by Phe 215, Gly290, Met388. 
The pocket is incompletely filled by the Rj'-iodo substiment. and accommodates the 
slightly larger 3'-isopropyI substiment by movement of the flexible Met388 side chain and 
the H7/H8 loop. This pocket can accommodate R3' substiments that are even larger than 

15 isopropyl, for example, a phenyl group. 

The above information will facilitate the design of high affmity agonists and 
antagonists by improving automated QSAR methodologies and informing manual 
modeling of pharmaceutical lead compounds. For example, the inclusion of discrete 
water molecules provides a complete description of hydrogen bonding in the polar pocket 

20 for use with pharmacophore development: also, the identification of mobile and immobile 
residues within the receptor suggests physically reasonable constraints for use in 
molecular mechanics/dynamics calculations. 

Example 12. Design of an Increased Affimty Ligand 
25 Direct interaction between the receptor and the ligand is limited in the polar 

pocket, which interacts with the Rj substiment. While the lack of complementarity may 
contain implications for biological regulation, it also provides an oppormnity for 
increasing affinity by optimizing the interaction between the amino acids of the polar 
pocket and the Ri substiment of a synthetic ligand. The strucmre of the receptor-ligand 
30 interactions described herein enables design of an increased affinity synthetic ligand 
having two complementary modifications: 
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1) Remove the positively charged amine. The strongly positive electrostatic 
potential predicted for the polar pocket suggests that the positively charged amine of the 
aminopropionic acid R, substituent may be detrimental to binding. Suitable groups for 
substitution are suggested by the nature of nearby hydrogen bond partners: for example. 
5 Thr 275 O or Ser 277 N. See e.g. Tables in Appendix 2. For example, any any 

negatively charged substituent would be compatible for interacting with the ammo acids of 
the polar pocket, including carboxylates, carbonyl, phosphonates. and sulfates, comprising 
0 to 4 carbons. Another example of an R, substitution is an oxamic acid that replaces the 
amine of the namrally occurring ligand with one or more carbonyl groups. 
10 2) Incorporate hydrogen bond acceptor and donor groups into the R,- 

substiment to provide broader interactions with the polar pocket scaffold. Such hydrogen 
bond acceptor and donor groups incorporated into the Rl -substituent will allow 
interactions that would otherwise occur with water molecules in the polar pocket. 
Specific waters include HOH 504 (hydrogen bonds with Ala 225 O and Arg 262 NH); 
15 and HOH 503 hydrogen bonds with Asn 179 ODl, Ala 180 N). both of which are present 
in all four complexes (TR LBD complexed with T3, TR LBC complexed with IpBr,. TR 
LED complexed with Dimit and TR LBD complexed with Triac). Analysis of the 
hydrogen bonding network in the polar pocket suggests replacement of HOH 504 with a 
hydrogen bond acceptor, and HOH 503 with an hydrogen bond donor (although the 
20 chemical namre of asparagine probably permits flexibility at this site). Thus, 

incorporating a hydrogen bond acceptor in an Rl substiment that could take the place of 
the HOH504 or incorporating a hydrogen bond acceptor in an Rl substituent that could 
positionally replace the HOH503, or a combination thereof, are methods of desigmng 

novel synthetic TR ligands. 
25 These two design approaches can be used separately or in combination to 

design synthetic ligands, including those in Table 4 (below). 

A corollary to this approach is to design specific interactions to the residues 
Arg262 and Asn 179. The goal is to build in interactions to these residues by desigmng 
ligands that have R, substiments that form hydrogen bonds with water molecules or 
30 charged residues in the polar pocket. 
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Table 4: Synthetic TR Ligands 




Rl R2 R3 R5 R6 X R'2 R*3 R'4 R*5 R'6 

10 C02H H Me Me H 0 H Me OH Me H 

CH2C02H II S Et SH El 

CH2CH2C02H Br Br nPr NH2 nPr 

CH2CH(NH2)C02H CI CI iPr iPr 

CK:H2C02H El El Pti riBu 

15 OCH2CH2C02H OH OH I "Pen 

NHCH2C02H NH2 NH2 Br nHcx 

NHCH2CH2C02H SH SH CI Ph 
CH2COCOC02H 

cycle 

20 NHC0C0C02H "T^' 
C0C02H 
CF2C02H 
COCH2C02H 



25 Any combination of ihe above substituenis in Uic biphcnyl ciher scaffold structure shown above may result in a 
potentially pharmacologically useful ligand for the thyroid hormone receptor. These novel ligands may be 
antagonists of the iliyroid receptor. 
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A strategy for designing synthetic ligands using the computational methods described 
herein is summed below: 



HO 




R20 



10 



For example. 




15 



20 



A = Hydrogen Bond Acceptor 

D = Hydrogen Bond Donor 

O = -OH, -CO 

RIO can be -OH, -CO 

R20 can be -CO 

R30 can be -COOH, -CONH2 

See also Table of synthetic TR Ligands 
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Dimit 

Data collection 

Cell dimensions 

a (A) 117.16 

b (A) 80 52 

c (A) 63.21 

P (•) 120.58 

Resolution (A) 2.2 

Obs. Reflections. 57031 
(no.) 

Unique Retleclions, 22327 
(no.) 

Completeness. (%) 87.0 



Phasing (15.0 - 2.5A) 

tR^, (%) 
No. of sites 
^Occupancy 
(Anomalous) 

§F„/E 

centric (acentric) 

15.0-5.0 A 

5.0-3.0 A 

3.0^2.5 A 
»Rc..J7c) 

15.0-5.0 A 

5.0-3.0 A 

3.0-2-5 A 

Mean figure of merit 0.62 

MR Phasing 
(10-3. 5A) 

Rotation Search: 

Eyler Angles (") 

§ correlation 
coefficient 

Translation Search: 

Fractional 
coordinate!; 
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TABLE 3: LBD- 122/4 10 

T3 IpBr, True 



1 17.19 

80.20 

63.23 

120.60 

2 0 

64424 

21023 

82.4 

3.5 



1 17.18 

80.12 

63.13 

120.69 

2.1 

66877 

23966 

93.7 
4.5 



1 18.19 

81 37 

63.73 

121.00 

2.45 

83573 

18453 

96.0 
7,5 



19.6 11.6 

3 2 

44.6 (19.8) 35.0 

50.2 (23.7) 35.0 
39.2 (22.3) 



3.67 (4.61) 
2.23 a.75) 
1.64 (1.99) 

33 
45 
60 



2.25 (3,09) 
1.25 (1.85) 
1.15 (1.57) 

44 

63 
65 



e, = 309.37 
Gj = 48.96 
Gj = 127.28 

34 3 

X = 0 1571 
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§ correlation 
Cocfficieni 

'R factor 

Refinenient 

* 

Resolution (A) 



1 Rtn« ( t> 

^ITt (t) 



/5.0'2.2 

20.5 
22.7 



5.0 - 2.0 

22.1 
24.0 



15.0 - 2.2 

21.4 
22.4 
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y = 0.000 
2 = 0.3421 

65.8 

31.2 
25-2.5 

23.6 
24.1 



10 



All publications and patent applications mentioned in this specification are herein 
incorporated by reference to the same extent as if each individual publication or patent 
application was specifically and individually indicated to be incorporated by reference. 
15 The nuclear receptor ligands, particularly the TR ligands. of these references are herein 
incorporated by reference and can be optionally excluded from the claimed compounds 
with a proviso. 

Headings and subheadings are presented only for the convenience of the reader 
and should not be used to construe the meaning of terms used within such headings and 
20 subheadings. 

The invention now being fully described, it will be apparent to one of ordinary 
skill in the art that many changes and modifications can be made thereto without 
departing from the spirit or scope of the appended claims. 
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Dimit 
Atom 

CIS 
CIS 
CIO 
C16 
CIS 
C19 

cia 

016 
CI9 
CS 

ce 

CIO 
C12 
01 
C13 

o-t:-: 

04 

C17 
03 

CIO 

C12 : 

CIO 

CM.: 

CIt 

C15 
N1 
CI 

cu 

C13 
CIS • 
C17 

oy 

04 

- C17 ^- 

04 
• N1 
C17 

Ni 

03 
" Hi ' 
NA 
Nt 
CS 

C19;; 

C7 
C9 

CII .* 

N1 
C13 

CIS 



Amino Acid 
In f ull length a 

215-PHE 

215- PHE 

216- PHE 
218-PHE 
218-PHE 

218- PHE 

219- THR 
22MLE 
222.ILE 
222-ILE 
222-tLE 
222-lLE 
222-ILE 
225-AUV 
22S-ALA 
228-ARG 
22a-ARG ■ ' 
22B-ARG 
22e-ARG:: 
25fl-MET 

■256-MeT.:"; ; 

256-MET 
256-MEX^ 
25».MeT 
259-MET: 
259-MET 
259-MET- 
2S9-MET 
259-MET: 
259-MET 
262-ARG: 
252- ARG 
262-ARG; " 

262- ARG 
: "262- ARG 

2«2.ARG 

: 2«3^ALA; 

263- ALA 

"* "iei ARCS 
:"2^ 

276^LEU 
;■ ■ 27W.EU . 
278-LEU 
27Q-LEU 
27fl-LEU 
278-LEU 
278-LEU 
27ft-LEU 
277.SER 
277-SER 



Ammo Acid 
Atom 

COl 
CE1 
0 
CB 

ca 

08 
CD2 
CG2 
CGI 
CD1 
GDI 
GDI 
GDI 
CDl 
CB 
CB 
CZ 
• NH2 
NH2 
NH2 
SO 

: SO 

CE 

■ CE 
0 

■ oi 
"o 

:''ol 

OS 
CB 
CB 

CB 
- CB 

; CB'' 

ca 

::-^-C0 ' 
CO 

"f- N:- 

CA 

CB ; 

NHl 

■''.O . 
CA 

c 

C01 

cor 

C02 
C02 

■ CD2 
C02 
H 



Distance 



N 



3,95 

3.33 

3.69 

3,89 

3.92 

4,13 

3.77 

3.63 

4.11 

4.18 

3.72 

3.53 

3.85 

4.13 

3.64 

3.90 
3.3Q 
3.58 
2.8S 
3.70 ' 

3.89;;^;;; 

3.88 

' -3-83^^; 
4.03 

3.42 

3.7t; 
4.20 

3-87,: Li. 
4.00 

4.01V ; V; 

3.5a__ 

■■3:i32:-r- 

3.85 

4.10^:-: 

3.61 

:;.3.7t:-.:" 
3.e9 
3.4a 

3.93 

3.62.::;-;.' 

3.51 

3.92-^:.. 

4 05 
'4.04 
4.09 
3.95 
4.13'^ 
4.17" 
4.14 
3.70 
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10 



Dimit 
Atom 



Amino Acia 
m full lengtn a 



Atom 



Distance A 



C17 

N1 

03 
C17 

03 
C13 

C7 
CIS 
C18 
C18 
C18 

C2 

C4 

C6 

C2 

C4 

02 

C4 

C6 

C8 

CIO 

01 

C20 

C8 

CIO 

01 

01 

ca 

CIO 
01 
C8 
CB 
01 
C15 
01 
C17 
03 
04 
C15 
C17 
N1 
03 
C17 
04 



277-SER 
277-SER 
277-S6R 
277-SER 
277-SER 
277-SER 
287-LEU 
290-GLY 

290- GLY 

291 - GUY 

292- LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
299-ILE 
381 -HIS 
381 -HIS 
381 -HIS 
381-HlS 
381 -His 
3B1-HIS 
381-HIS 
388-MET 
401-PHE 
401 -PHE 
401-PHE 
401-PHE 

3-H3O 
502 

3- HiO 

2-HiO 
2-HjO 
503 

4- HiO 
504 



N 
N 
N 

CA 
CA 
OG 
C02 
C 

o 

CA 
N 

CG 
CG 
CG 
CD! 
CDl 
GDI 
CD2 
CD2 
CD2 
CD2 
CD2 
CDl 
CD2 
CD2 
G02 
CE1 
NE2 
NE2 
NE2 
CE 
CEl 
CEl 
CZ 
CZ 
01 
01 
01 
01 
01 
01 
01 
01 
01 



3.69 

3.30 

3.19 

3.92 

3.35 

3.92 

3.90 

4.04 

3.54 

4.04 

4.20 

4.18 

3.88 

4.01 

3.88 

4.02 

4.07 

4.05 

3.72 

3.69 

3-98 

4.16 

3.87 

3.90 

3.84 

3.40 

3.72 

3.47 

3.51 

2- 64 
3.90 
4.19 
3.37 

3- 97 
3.28 
3.77 
3.13 
3.72 
4.04 
3.84 
3.35 
2.56 
3.92 
2.72 



Legend to Table XX. The Ubie lists the interactions with Dimit (DMT), The column headings 
50 are as follows: 

The atom of Dimit that interacts with The ammo acid of the receptor. These are also 
numbered in figure X. 

U2 The ammo acid in the full length rTRa that interacts with the hgand. 
ifZ The name of the atom in the ammo acid Istandard nomenclature! where the interaction 
55 occurs. 

#4 The distance m A between Dimit and the orotem atom. 



J 



APPENDIX : 



wo 97/21993 



78. 



PCT/US96/20778 



Triac 


Amino Acid 


Amino Acid 




Atom 


in full length a 


Atom 




W 1 


2 1 8-PHE 


0 


TCI 


11 


221 -ILE 


CDl 


/I 1 C 

* . 1 D 


11 


222-ILE 


CA 


A 1 c 
4. 1 D 


n 


222-ILE 


CB 


4.U J 


11 


221 -ILE 


CG 1 


3.92 


C8 


222-ILE 


CDl 


4. 1 2 


CIO 


222-ILE 


CD! 


3,77 


C12 


222-ILE 


CDl 


3.79 


C13 


225-ALA 


CB 


4.17 


C3 


225*ALA 


CB 


3.86 


CIO 


25a-MET 


SD 


3.45 


012 


25a-MET 


SD 


3.73 


CIO 


258-MET 


CE 


3.66 


C12 


256-MET 


CE 


3.77 


Cll 


256-MET 


C 


3.68 


Cil 


256-MET 


0 


3.24 


CI 


259-MET 


0 


3.93 


C11 


259-MET 


0 


3.24 


03 


259-MET 


0 


4.09 


CI 


259-MET 


CB 


3.89 


C13 


259-MET 


0 


3.74 


C14 


259-MeT 


0 


3.96 


CI 


259-MET 


CB 


3.89 


Cll 


259-MET 


CB 


3.68 


C13 


259-MET 


CB 


4.01 


C14 


262-ARG 


CB 


4.07 


04 


262-ARG 


CB 


3.60 


C17 


282-ARG 


CD 




04 


262-ARG 


CD 




03 


263-ALA 


N 


3.79 


03 


263-ALA 


CA 


3.67 


03 


283-ALA 


CB 


3.49 


cn 


107-ALA 


CB 


4.00 


03 


107-ALA 


CB 


3.49 


03 


266-ARG 


NHl 


3.00 


03 


275-THR 


0 


3.20 


03 


278-LEU 


CA 


3.1 1 


N1 


276-LEU 


C 


3.52 


03 


1 20-LEU 


N 


4.04 


C14 


1 20-LEU 


CA 


3.98 


03 


120-LEU 


CA 


3.11 


C14 


1 20-LEU 


C 


3.98 


03 


120-LEU 


CB 


3.95 


02 


276-LEU 


CD! 


4.03 


11 


276-LEU 


CDl 


4.10 


C7 


276-LEU 


CD2 


3.84 


C9 


276-LEU 


CD2 


3.73 


CM 


276-LEU 


CD2 


4.06 


02 


276-LEU 


CD2 


4.10 
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15 



20 



25 



30 



( nac 
Atom 



C13 

C14 

04 

03 

C14 

04 
C3 
12 
12 
C4 
C6 
C2 
C4 
12 
C4 
C6 
C8 
CIO 
01 
13 
CB 
CIO 
01 
01 
C8 
CIO 
01 
01 
01 
01 



Ammo AciO 
in full lengtn a 
277-SER 
277-SER 
277-SER 
277-SER 
277-SER 
277-SER 
277-SER 
29Q-GLY 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
292-LEU 
299-ILE 
381-HlS 
381-HlS 
381-HIS 
381-HlS 
381-HIS 
381-HIS 
381-HIS 
388-MET 
401-PHE 
401 -PHE 



Amino Acid 
Arom 



Distance A 



lM 

N 
N 
N 

CA 
CA 
OG 
□ 

CG 
CG 
CG 
CDl 
CD1 
CD2 
CD2 
CD2 
CD2 
CD2 
CD2 
CDl 
CD2 
CD2 
G02 
CEl 
NE2 
NE2 
NE2 
CE 
CEl 
CZ 



4.06 

3.13 

3.28 

3.05 

3.76 

3.52 

3.87 

3.57 

3-94 

3.95 

3.65 

4.1 1 

3.85 

3.98 

4.1 1 

3.44 

3.28 

3.88 

3.35 

3.77 

3.87 

3.90 

3.20 

3.81 

3.51 

3.52 

2.64 

4.03 

3.86 

3.70 



Legend to Table XX. The table lists the .nteract.ons with tr.ac. The column headmgs are as 



follows: 



U^ The atom of tr.ac that .nteracts w.th the amino ac.d of the receptor. These are also 

35 numbered in figure X. , 

tt2 The amino acid in the full length rTRo that interacts w.th the l.gand. 

A The name of the atom m the amino ac.d (standard nomenclature) where the .nteract.on 



occurs. 

#4 The distance in A between triac and the protein atom. 
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f_= — _ 

ipBR; Atom 


Amino Actd 
lO full length a 


Amino Acid 
Atom 


Distance 


C16 


2 1 5-PHE 


GDI 


A ^1 

4 U 1 


CI 6 


2 1 5-PHE 


CEl 


3, 78 


BR1 


218-PHE 


0 


3 24 


BRl 


218-PHE 


c 


3,98 


C16 


218-PHE 


CB 


3,81 


C18 


218-PHE 


CB 


3.92 


BRl 


218-PHE 


CB 


4 08 


CIS 


218-PHE 


CD2 


3 92 


CIS 


219-THR 


CG2 


3 45 


BRl 


22MLE 


CGI 


3 81 


BRl 


221-ILE 


CD1 


4 07 


C6 


222-tLE 


CD1 


4 07 


C8 


222-ILE 


CD1 


3 64 


CIO 


222-ILE 


CD1 


3.50 


C12 


222-ILE 


CD1 


3.82 


01 


222-ILE 


GDI 


4 08 


C13 


225*ALA 


CB 


3.76 


04 


225-ALA 


CB 


4.01 


04 


22B-ARG 


cz 


3.92 


C17 


22B-ARG 


NH2 


3.26 


03 


228-ARG 


NH2 


3.43 


04 


228-ARG 


NH2 


2.79 


CIO 


256-MET 


SO 


3.55 


C12 


256-MET 


SD 


3-71 


CIO 


256-MET 


CE 


3.90 


C12 


256-MET 


CE 


3.75 


BR2 


256-MET 


CE 


4.03 


CI 1 


259-MET 


L 




CI 1 


259-MET 


0 


3.52 


C15 


259-MET 


0 


3.44 


Nl 


259-MET 


0 


3.76 


C11 


259-MET 


CB 


3.87 


C15 


262-ARG 


CB 


4 03 
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Ir 


IpBRj Atorri 


Amino Acid 
m full tength c 


Amino Actd \ 
Atom i 


Distance 
A 


r 


r 1 7 


262-ARG 


1 

CB 1 


3.56 






262*ARG 


CB 


3.55 




04. 


262-ARG 


CB 


3.91 




n 7 


262-ARG 


CD 


4.09 


c 

J 


04. 


262-ARG 


CD 


3.71 




hi 1 
IM 1 


263-ALA 


N 


3.61 




(N 1 


263-ALA 


CA 


3.59 




N 1 


263-ALA 


CB 


3.54 




U J 


268-ARG 


NHl 


3.93 


1 f\ 


M 1 


275-THR 


0 


3.43 




N 1 


276-LEU 


CA 


3.46 




W 1 


276-LEU 


C 


3.83 




Cb 


276-LEU 


CD! 


4.02 






276-LEU 


CD2 


4.00 


1 J 




276-LEU 


CD2 


3.81 




Li 1 


276-LEU 


CD2 


3.91 




U 1 O 


277-SER 


N 


3.79 




C 1 b 


277-SER 


N 


3.63 




Ci 7 


277-SER 


N 


3.70 


2\J 


N 1 


277-SER 


N 


3.17 




U J 


277-SER 


N 


3.37 




/"IT 


277-SER 


CA 


3.89 




U J 


277-SER 


CA 


3.43 




L. 1 o 


277-SER 


OG 


3.66 






287-LEU 


CD! 


4.05 




r" 1 ft 


290-GLY 


C 


4.04 




1 ft 


290-GLY 


0 


3.48 




1 ft 


291 -GLY 


CA 


4.02 




C4 


292-LEU 


CG 


3.89 


30 


C6 


292-LEU 


CG 


4.02 




C2 


292-LEU 


CD! 


3-79 




C4 


292-LEU 


CD! 


3.96 




02 


292-LEU 


CD! 


3.97 
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IdBR. Atom 


Ammo Acid 
in full length a 


Amino AciC 
Atom 


Distance 
A 






292-LEU 




4.07 


1 






292-LEU 


rn 9 


3 75 




Co 




292-LEU 


L.L/Z 


3 67 




CIO 




292-LEU 




3 92 




BR2 




299-ILE 


CU I 


3 68 




C8 




361 -HIS 

O V 1 1 1 1 W 


^ 1 

CD2 


3 67 




CIO 




381 -HIS 


CD2 






01 






GD2 


3 




01 




JO 1 * n I o 


CE1 


J D ^ 




C8 




OQ1 .UlC 

JO 1 *nFo 


NE2 


J. JO 




CIO 






NE2 


0 OH 




01 






NE2 






C8 




401 -PHF 


CE1 


4 07 




01 






CE1 


3 19 




C16 






CZ 






01 




Am PUT 


cz 






C17 






01 






03 




502H,O 


01 


3.40 




04 




3-HjO 


01 






f 1 c 

L. 1 3 




2-HjO 


m 






C17 




2-H2O 


01 






Nl 




502HjO 


01 


3.12 




03 




503H3O 


01 


2.27 




017 




4-H2O 


01 






04 




504HjO 


01 


1.78 





Legend to Table XX. The table lists the interactions with lpBr2. The column headings are as 
foltows: 

#1 The atom of lpBr2 that interacts with the amino acid of the receptor. These are also 
30 numbered in figure X. 

U2 The ammo acid in the full length rTRo that interacts with the ligand. 

#3 The name of the atom in the amino acid (standard nomenclature! where the interaction 

occurs. 

^4 The distance in A between lpBr2 and the protein atom. 
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10 



15 



20 



25 



30 



'3 Atom 



12 

n 
n 

C4 

n 
11 

C8 
CIO 
C12 
01 
C13 
C3 
04 
04 
C15 
03 
04 
CIO 
C12 
C10 
C12 

cn 
cn 

C14 
N1 

ci 

C11 
C13 
C15 
C14 
03 
04 
04 



Amino Acid 
in full length a 



215-PHE 
218-PHE 
21 8-PHE 
218-PHE 
218-PHE 

221 - ILE 

222- ILE 
222-ILE 
222-tLE 
222-ILE 
225-ALA 
225-ALA 
228-ARG 
228-ARG 
228-ARG 
228-ARG 
228-ARG 
256-MET 
256-MET 
256-MET 
256-MET 
259-MET 
259-MET 
259-MET 
259-MET 
259-MET 
259-MET 
259-MET 
262-ARG 
262-ARG 
262-ARG 
262-ARG 
282-AHG 



Amino Acid 
Atom 



GDI 
0 

CB 
CB 
CB 
CGI 
CD1 
CD! 
CD! 
CD1 
CB 
CB 
NHl 
CZ 
NH2 
NH2 
NH2 
SD 
SD 
CE 
CE 
C 
0 
O 
0 

CB 
CB 
CB 
CB 
CB 
CB 
CB 
CD 



Distance 



4.08 

3.19 

3.99 

4.04 

3.99 

4.01 

3.99 

3.99 

3.57 

3.66 

3.66 

4.04 

3.23 

3.45 

3.54 

3.90 

2.86 

3.73 

3.90 

3.97 

3.92 

3.95 

3.59 

3.51 

3.88 

4.06 

3.77 

3.96 

3.61 

4.02 

3.65 

3.92 

3.72 
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j T3 Atom 


Ammo Acid 

in f 1 1 II IpnnT h rr 


Amfno Acid 
Atom 


I 

Distance 
A 


M 1 

IN 1 


T Al a 
D J-MLM 




3.81 


M 1 




PA 
CM 




N 1 


T CO A 1 A 


PR 


3 63 


N 1 


2 7b- 1 HH 


L/ 
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n T o r c 1 1 
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P A 
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N 1 
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276-LEU 


CD! 


4 , Ub 


C7 
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3. HO 
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277-SER 
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277-SER 
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NHl 


3.23 
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277-SER 


N 


J.U / 
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N 
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CA 
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O. / / 
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277-SER 


CA 


J. J4 
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277-SER 


OG 




12 


290-GLY 




J. 3U 


C4 


292-LEU 


L.Li 


3.95 


C8 


292-LEU 




1 Ul 
O.O J 


C2 


292-LEU 


CU 1 




C4 


292-LEU 


p m 




C4 


292-LcU 


pn 9 




C6 


2y2-LtU 






C8 


292-LcU 


pn "5 




L. 1 u 


^ 3 ^ * l.C^ w 


CD2 


3.96 


01 


292-LEU 


CD2 


3.71 


13 


299-ILE 


CDl 


3.74 


C8 


381-HIS 


CD2 


3.94 


CIO 


381-HIS 


CD2 


3.97 


01 


381-HIS 


CD2 


3.39 
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T3 Atom 


Ammo Acid 
m full length a 


Amino Actd 
Atom 


Distance 
A 





01 


3B1-HIS 


CDl 


3.82 






C8 


381-HIS 


NE2 


3.47 






CIO 


381-HlS 


NE2 


3.55 






01 


381-HIS 


NE2 


2.70 
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388-MET 
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3.88 






01 


401 -PHE 
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3.52 






01 


401 'PHE 
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3.32 






03 


502 


01 


2.51 






04 


3-HiO 


01 








N1 


2-H3O 


01 








03 


503 


01 


2.81 






04 


504 


01 


2.73 





10 



Legend to Table XX. The table lists the interactions with T3. The column headings are as 
15 follows: 

#1 The atom of T3 that interacts with the amino acid of the receptor. These are also 
numbered in figure X. 

#2 The amino acid in the full length rTRa that interacts with the ligand. 

#3 The name of the atom in the amino acid (standard nomenclature! where the interaction 

20 occurs. 

#4. The distance in A between T3 and the protein atom. 
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Coordination S;njcrjre of Thyroid Hormone Receptor and T3 
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90. 

WHAT IS CLAIMED IS: 

1. A crystal of a TR LBD comprising: 

1) a TR LBD and a TR LBD ligand or 

2) a TR LBD without a TR LBD ligand; 

5 wherein said crystal diffracts with at least a 2.0 to 3.0 angstrom resolution and has 

a crystal stability within 5% of its unit cell dimensions. 

2. The crystal of claim 1, wherem said TR LBD has at least 200 amino acids. 

10 3, The crystal of claim 2, wherein said TR LBD is from a TR protein selected from 
the group consisting of TR amino acid sequence 122 to 410 of rat TR-a, 157 to 410 rat 
TR-a and 202 to 461 of human TR~{3. 

4. The crystal of claim 2, wherein said TR LBD ligand is a compound being of a 
15 formula: 



20 




25 



wherein Ri is 

-O-CH.CO.H. -NHCH.COjH, 
30 -COM. -CH^CO.H, -CH.CH.CO.H. -CH.CHXH^CO.H, 





PCT/L'S96/20778 
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-CH3CH(NH,)CO,H. -CH,CH[NHCOCH<^)JCO>H. -CH,CH[NHCO(CH,),sCH, 
]CO:H. -CH,CH[NH-FMOC]CO:H, -CH.CHlNH-tBOClCO^H, or a carboxylate 
connected to the ring with a 0 to 3 carbon linker. 



-PO3H,, -CH^POjH,, -CH,CH,P03H,. -CH,CHNH,P03H,, 
-CHjCH[NHCOCH<^,]P03H„ -CH:CH[NHCO(CH2),3CH3]P03H,, 
-CH,CH[NH-FM0C]P03H,, -CH, CH[NH-tB0ClP03H;, or a phosphate or 
phosphonate connected to the ring with a 0 to 3 carbon linker, 

-SO3H, .CH2SO3H. -CHXH,S03H, -CH.CHNHjSOjH. 
-CH,CH[NHC0CH(i>JS03H, -CH,CH[NHCO(CH,),5CH3lS03H. 
-CH:CH[NH-FMOC]S03H, CU. CH[NH-tBOClS03H. or a sulfate or sulfite 
connected to the ring with a 0 to 3 carbon linker, 



or acts as the functional equivalent of CH2CH(NH,)C0iH of T3 in the molecular 
recognition domain when bound to a TR, wherein said R, can be optionally 
substituted with an amine. 



wherein R, is 

20 

halogen, CF3, OH. NH„ SH, CH3, -Et, 
or acts as the functional equivalent of H in the molecular recognition domain when 

bound to a TR, 



25 wherein R3 is 



-H, halogen, -CF,, -OH, -NH^, -N,. -SH, -CH3, -Et. 

or acts as the functional equivalent of I in the molecular recognition domain when 
bound to a TR, 

wherein R5 is 
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-H, halogen, -CF,. -OH. -NH., -N3, -SH. -CH3, -Et. or acts as the functional 
equivalent of I in the molecular recognition domain when bound to a TR, and Ri 
can be identical to R., 



5 wherein R^ is 



-H, halogen. -CF3, -OH, -NH,, -SH, -CH3, or acts as the functional equivalent of 
H in the molecular recognition domain when bound to a TR» and R. can be 
identical to R^, 

10 

wherein R2' is 

-H, halogen, -CF3. -OH. -NHj. -N3, -SH, -CH3, -Et, or acts as the functional 
equivalent of H in the molecular recognition domain when bound to a TR. 

15 

wherein R,' is any hydrophobic group, including 

halogen, -CF3, -SH. alkyl, aryl, 5- or 6-membered heterocyclie, cyano, or acts as 
the functional equivalent of I in the molecular recognition domain when bound to a 
20 TR, 

wherein R4' is 



-H, halogen, -CF3, -OH. -NHj. NH3, -N(CH3)3, carboxylate. phosphonate, 
25 phosphate or sulfate. -SH. -CHj, -Et, or akyl, aryl or 5- or 6-membered 

heterocyclic aromatic attached through urea or carbamate linkages to O or N or S 
at the R,' position, or acts as the functional equivalent of OH in the molecular 
recognition domain when bound to a TR, 



30 wherein R/ is 





10 
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-H. -OH, -NH,, -N(CH3): -SH -NH3, -N(CH3)3. carboxylate. phosphonate. 
phosphate, sulfate, branched or straight chain alkyl having 1 to 9 carbons, 
substituted or unsubstituted aryl, wherein said substituted aryl is substituted with 
halogen or I to 5 carbon alkyl and wherein said aryl is optionally connected to the 
ring by a -CH^-, aromatic heterocycle having 5 to 6 atoms, wherein said 
heterocycle may be substituted with one or more groups selected from -OH. -NH„ 
-SH, -NHj, -N(CH3)3, carboxylate, phosphonate, phosphate, sulfate, heteroalkyU 
heteroaryl, arylalkyl, heteroaryl alkyl, polyaromatic, polyheteroaromatic, wherein 
said R5' may be substituted with polar or charged groups, 

wherein R^' is 



-H. halogen. -CF3, -OH. -NH^, -SH, -CH3, -Et, or acts as the functional 
equivalent of H in the molecular recognition domain when bound to a TR. 

15 

wherein X is 

O, S, SO:, NH, NR„ CHj, CHR,. CR,R,. wherein R, is alkyl, aryl or 5- or 
6-membered heterocyclic aromatic, 
20 and wherein said TR LBD ligand has an apparent Kd for binding TR LBD of 1 /xM or 

less. 

5 The crystal of claim 2, wherein said crystal of TR LBD ligand has the following 
unit cell dimensions m angstroms: a= 117.00 ± 2%. b = 80.00 ± 2%, c = 63.00 ± 2%. 
25 at a /3 angle of 120.00 ± 2% degrees and a monoclinic space group C2. 

6. The crystal of claim 5, wherein said crystal coordinates said TR LBD ligand and 
said crystal is selected from the group consisting of crystals having one of the following 
sets of properties: 

30 
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1) unit cell dimensions in angstroms: a = 1 17. 16, b = 80.52, c =63.21. at a 
angie of 120.58 degrees, a monoclinic space group C2, 2.2 angstrom resolution 
and completeness value of 87.0%, 

5 2) unit cell dimensions in angstroms; a=117.19, b = 80.20, c =63.23, at a /? 

angle of 120.60 degrees, a monoclinic space group C2, 2.0 angstrom resolution 
and completeness value of 82.4%, and 

3) unit cell dimensions in angstroms; a = 1 17.18, b = 80.12, c = 63.13, at a P 
10 angle of 120.69 degrees, a monoclinic space group C2, 2.2 angstrom resolution 

and completeness value of 93.7%. 

7- The crystal of claim 2, wherein said crystal further comprises a human protein. 

15 8. The crystal of claim 7, wherein said crystal has the following unit cell dimensions 
in angstroms: a = b = 68.448 ± 2% and c= 130.559± 2% at an a angle of 90.00", a (S 
angle of 90.00° and a j angle of 120.00° and having a trigonal space group p3(l)21. 

9. The crystal of claim 8, wherein said crystal further comprises a human protein. 
20 10. A computational method of designing a nuclear receptor synthetic ligand 
comprising: 

1) using a three dimensional model of a crystallized protein comprising a nuclear 
receptor LBD with a bound ligand to determine at least one interacting amino acid of the 
nuclear receptor LBD that interacts with at least one first chemical moiety of said bound 

25 ligand, and 

2) selecting at least one chemical modification of said first chemical moiety to 
produce a second chemical moiety with a structure that either decreases or increases an 
interaction between said interacting amino acid and said second chemical moiety 
compared to said interaction between said interacting amino acid and said first chemical 

30 moiety. 



11 



The method of claim 10. wherein steps 1 and 2 are repeated. 





PCT/US96/20778 

WO 97/21993 

95. 

12. The method of claim 10, further comprising generatmg said three dimensional 
model of said crystallized protein comprismg a TR LBD with a bound TR ligand. 

13. The method of claim 12, wherein said three dimensional model is generated by 
5 comparing isomorphous ligand derivatives to produce improved phasing. 

14. The method of claim 13, wherein said isomorphous ligand derivatives are provided 
by substituting at least one of the R5, R3, 'R5 and 'R3 positions of a thyronine derivative 
with Br or I. 

10 

15. The method of claim 14, further comprising determining a change in interaction 
between said interacting ammo acid and said ligand after chemical modification of said 
first chemical moiety. 

15 16- The method of claim 12, wherein said selecting uses said first chemical moiety 
that interacts with at least one of the said interacting amino acids listed in APPENDIX 2. 

17. The method of claim 16, wherein said chemical modification enhances hydrogen 
bonding interaction, charge interaction, hydrophobic interaction, Van Der Waals 

20 interaction or dipole interaction between said second chemical moiety and said interacting 
amino acid compared to said first chemical moiety and said interacting amino acid. 

18. The method of claim 17, wherein said chemical modification reduces hydrogen 
bonding interaction, charge interaction, hydrophobic interaction. Van Der Waals 

25 interaction or dipole interaction between said second chemical moiety and said interacting 
amino acid compared to said first chemical moiety and said interacting amino acid. 

19. The method of claim 18, wherein said first chemical moiety is at least 4.5 
angstroms away from at least one distant amino acid and said distant amino acid is not 

30 any of the said interacting amino acids listed in claim 17. 
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20. The method of claim 18, wherein said first chemical moiety is 6 to 12 angstroms 
away from a distant amino acid. 

21. The method of claim 19, wherein said chemical modification extends toward said 
5 distant amino acid and produces hydrogen bonding interaction, charge interaction, 

hydrophobic interaction, Van Der Waals interaction or dipole interaction between said 
second chemical moiety and said distant amino acid. 

22. The method of claim 19, wherein said chemical modification sterically hinders said 
10 distant amino acids from binding other ammo acids but does not sterically hinder said 

interacting amino acids. 

23. The method of claim 17 or 19 further comprising determining a change in 
interaction between said interacting amino acid and said hgand after said chemical 

15 modification, wherein said determining and modifying are performed using a computer 
program to represent chemical structures of said interacting amino acid and ligand. 

24. The method of claim 22, wherein said chemical modification is at a R5' position of 
a thyronine derivative. 



25. The method of claim 24, wherein said chemical modification sterically hinders 
activation helix function. 

26. The method of claim 25, wherein said chemical modification does not interfere 
25 with an interaction between said interacting amino acid and an atom from R 1 to 6 

positions of said ligand. 

27. The method of claim 26, wherein said activation helix is helix H12. 



20 



30 28. The method of claim 27. wherein said chemical modification comprises a planar 
structure, equivalent in length to at least a 3 carbon alkyl that projects out from the plane 
of the prime ring of said thyronine derivative at least 30"*. 
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29. The method of claim 10 further comprismg generating said three dimensional 
model of a nuclear receptor other than TR using said nuclear receptor's LBD amino acid 
sequence and using said crystallized protein comprising a TR LBD with a bound TR LBD 
ligand. 

5 

30. The method of claim 29, wherein said chemical modification enhances hydrogen 
bonding interaction, charge interaction, hydrophobic interaction. Van Der Waals 
interaction or dipole interaction between said second chemical moiety and said interacting 
amino acid compared to said first chemical moiety and said interacting amino acid. 

10 

31. The method of claim 30, wherein said chemical modification reduces hydrogen 
bondmg interaction, charge interaction, hydrophobic interaction. Van Der Waals 
interaction or dipole interaction between said second chemical moiety and said interacting 
amino acid compared to said first chemical moiety and said interacting amino acid. 

15 

32. The method of claim 31, wherein said chemical modification is only of said first 
chemical moiety that is at least 3 angstroms away from at least one distant amino acid and 
said distant amino acid is different from said interacting amino acid. 

20 33. The method of claim 33, wherein said first chemical moiety is 4 to 12 angstroms 
away from a distant amino acid. 

34. The method of claim 33. wherein said chemical modification extends toward said 
distant amino acid and produces hydrogen bonding interaction, charge interaction, 

25 hydrophobic interaction, Van Der Waals interaction or dipole interaction between said 
second chemical moiety and said distant amino acid. 

35. The method of claim 34, wherein said chemical modification sterically hinders said 
distant amino acids from binding other amino acids but does not sterically hinder said 

30 interacting amino acids. 
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36. The method of claim 35, wherein said first chemical moiety comprises C-H or 
C-OH- 

37. The method of claim 36. wherein said second chemical moiety comprises a longest 
entity equivalent to at least 4 double bonded carbons in length and a volume of at least 4 

5 double bonded carbons. 



38. A computational method of designing a nuclear receptor antagonist from a nuclear 

receptor agonist comprising: 

1) determining a structure of a molecular recognition domain of said agonist using 
10 a three dimensional model of a crystallized protein comprising a nuclear receptor LBD, 

and 

2) selecting at least one chemical modification of said agonist that provides a 
ligand structure that extends beyond a binding site for said agonist and in ihe direction of 
at least one protein domain important in nuclear receptor function. 

15 

39. The method of claim 38, wherein said protein domain is either: 

a) a transcription activation domain of said LBD, 

b) a repressor binding domain of said LBD, 

c) a DNA binding domain of said nuclear receptor, 

20 d) a heat shock protein binding domain of said nuclear receptor. 

e) a dimerization domain of said LBD, or 
0 a hinge region to said DNA binding domain. 

40. The method of claim 38, wherein said crystallized protein comprises said nuclear 
25 receptor LBD bound to a nuclear receptor ligand. 

41. The method of claim 40 further comprising generating said three dimensional 
model of said crystallized protein comprising a TR LBD with a bound TR ligand. 



30 42. The method of claim 40 wherein said three dimensional model is for a TR LBD. 
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43. The method of claim 41 farthei comprising generating said three dimensional 
model of a nuclear receptor other than TR using said nuclear receptor's LBD amino acid 
sequence. 

5 44. The method of claim 39 wherein said three dimensional model is for an ER LBD. 

45. The method of claim 42 wherein said LBD is from a receptor selected from a 
group consisting of glucocorticoid receptor, estrogen receptor, retinoid receptor and 
vitamin D receptor. 

10 

46. The method of claim 39, wherein said chemical modification mimmizes a loss of 
hydrogen bonding interaction, charge interaction, hydrophobic interaction. Van Der 
Waals interaction or dipole interaction between said molecular recognition domain and 
said antagonist compared agonist hydrogen bonding interaction, charge interaction, 

15 hydrophobic interaction or dipole interaction with said molecular recognition domain. 

47. The method of claim 39, wherein said chemical modification reduces hydrogen 
bonding interaction, charge interaction, hydrophobic interaction or dipole interaction 
between said molecular recognition domain and said antagonist compared with the agonist 

20 hydrogen bonding interaction, charge interaction, hydrophobic interaction. Van Der 
Waals interaction or dipole interaction with said molecular recognition domain but said 
chemical modification still permits said antagonist to bind to said nuclear receptor with an 
affinity of lOOnM Kd or less. 

25 48. The method of claim 39 or 47, wherein said molecular recognition domain is at 
least 3 angstroms away from at least one distant amino acid and said distant amino acid 
does not contribute significantly to binding said molecular recognition domain. 



49. The method of claim 48, wherein said first molecular recognition domain is 4 to 
30 12 angstroms away from a distant amino acid. 
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50. The method of claim 49. wherein said chemical modification extends toward said 
distant amino acid and produces a hydrogen bonding interaction, charge interaction, 
hydrophobic interaction, Van Der Waals interaction or dipole interaction between said 
distant amino acid and said molecular recognition domain. 

5 

51. The method of claim 50, wherein said chemical modification sterically hinders said 
distant amino acids from binding other amino acids but does not sterically hinder amino 
acids that bind said molecular recognition domain. 



10 52. The method of claim 51, wherein said chemical modification of said ligand is at a 
C-H or C-OH. 

53. The method of claim 52, wherein said second chemical modification consists of a 
longest entity equivalent to at least 4 double bonded carbons in length and has volume of 

15 at least 4 double bonded carbons in length. 

54. The method of claim 44, wherein said three dimensional model is represented in 
two dimensions- 

20 55. A computational method of designing a nuclear receptor super agonist or 

antagonist comprising: 

1) determining at least one interacting amino acid of a nuclear receptor LBD that 
interacts with at least one first chemical moiety of said ligand using a three dimensional 
model of a crystallized protein comprising a nuclear receptor LBD with a bound ligand, 

25 and 

2) selecting at least one chemical modification of said first chemical moiety to 
produce a second chemical moiety with a structure to reduce or enhance an interaction 
between said interacting amino acid and said second chemical moiety compared to said 
interaction between said interacting amino acid and said first chemical moiety. 

30 

56. The method of claim 55, wherein said chemical modification enhances hydrogen 
bonding interaction, electrostatic interaction, charge interaction, hydrophobic interaction. 
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Van Der Waals interaction or dipoie interaction between said second chemical moieiy and 
said interacting amino acid compared to said first chemical moiety and said mteraciing 
amino acid. 

5 57. The method of claim 55, wherein said chemical modification changes a 

carboxylate moiety of said first chemical moieiy to a phosphonate or phosphate to make 
said second chemical moiety . 

58. The method of claim 55, wherein said nuclear receptor is TR and said chemical 
10 modification enhances said interaction between said second chemical moiety and at least 
one of the following arginines: Arg 262, Arg 266 or Arg 228 of the rat a-TR or an 
arginine of human a-TR or 0-JR that corresponds in its three dimensional position in said 
three dimensional model to either said arginines: Arg 262, Aig 266 or Arg 228. 

15 59. The method of claim 58. wherein said chemical modification changes a side group 
to nil space normally occupied by water when an agonist or naturally occurring ligand is 
bound to said LBD. 

60. The method of claim 59, wherein said chemical modification snugly fits into a 
20 pocket or surface of a binding cavity of said LBD and complements to a charge or 

hydrophobicity or both of said pocket or surface. 

61. The method of claim 60, wherein said chemical modification changes said first 
chemical moiety either a first carbon bonded either to 1) a first phenyl at Rl of T3 or 2) 

25 a first phenyl at Rl of a T3 agonist with a T3 carboxylate group to a second chemical 
moiety a two carbon space group. 

62. The method of claim 55, wherein said chemical modification changes said first 
chemical moiety either 1) a carbon linking a first phenyl at R4 of T3 to a second phenyl 

30 ring of T3 or 2) a carbon linking a first phenyl at R4 of a T3 agonist with a second 
phenyl ring of said T3 agonist to a second chemical moiety mono- or geminal-substiruted 
carbon group. 
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63. The method of claim 55, wherein said chemical modification changes said first 
chemical moiety either 1) a carbon linking a first phenyl at R4 of T3 co a second phenyl 
ring of T3 or 2) a carbon linking a first phenyl at R4 of a T3 agonist with a second 
phenyl ring of said T3 agonist. 
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64. A compound of the formula: 




wherein Rj is 
15 -O-CH2CO3H. -NHCH2CO2H, 

^CO,H, -CH.CO^H, -CH^CH^CO^H, -CH,CHXH,CO,H, 
-CH,CH(NH,)C02H, -CHXH[NHCOCH0JCO,H, -CH,CH[NHCO(CH,),5CH3 
ICOjH, -CHXH[NH-FM0C1C03H, -CH,CH[NH-tBOC]COjH. or a carboxylate 
connected to the ring with a 0 to 3 carbon linker, 

20 

-PO3H2, -CH,P03H„ -CHXHjPOjHj, -CH,CHNH,P03H,. 
-CH:CH[NHCOCH<^2]P03H,. -CH,CH[NHCO(CH,),5CH3lP03H,, 
-CH,CH[NH-FMOC]P03Hj, -CH, CHINH-tBOCJPOjH,, or a phosphate or 
phosphonate connected to the ring with a 0 to 3 carbon linker, 

25 

-SO3H. -CH2SO3H, -CH.CH^SOjH, -CHXHNH.SOjH. 
-CHiCH[NHC0CH<^),]S03H. -CH2CH[NHCO(CH2),5CH3]S03H, 
-CH,CH[NH-FM0C1S03H. -CH, CH[NH-tB0C]S03H, or a sulfate or sulfite 
connected to the ring with a 0 to 3 carbon linker, 
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or acts as the functional equivalent of CH.CHCNHOCO.H of T3 in the molecular 
recognition domain when bound to a TR, wherein said R, can be optionally 
substituted with an amine, 

5 wherein R, is 

H, halogen, CF3. OH, NH^, SH, CH3. -Et, 

or acts as the functional equivalent of H in the molecular recognition domain when 
bound to a TR, 

10 

wherein R3 is 

-H. halogen, -CF3, -OH, -NH^, -N3, -SH. -CH^, -Et, 

or acts as the functional equivalent of I in the molecular recognition domain when 
15 bound to a TR, 

wherein R5 is 

-H, halogen, -CF3. -OH, -NH„ -N3, -SH, -CH3, -Et, or acts as the functional 
20 equivalent of I in the molecular recognition domain when bound to a TR, and R3 

can be identical to R5, 

wherein Kf, is 

25 -H, halogen, -CF3, -OH, -NH,, -SH, -CH3, or acts as the functional equivalent of 

H in the molecular recognition domain when bound to a TR, and R. can be 
identical to R^j, 

wherein R,' is 

30 

-H, halogen. -CF3, -OH, -NH^, -N3. -SH, -CH3. -Et, or acts as the functional 
equivalent of H in the molecular recognition domain when bound to a TR, 
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wherein Rj' is any hydrophobic group, including 

halogen, -CF,, -SH, alkyl, aryl, 5- or 6-membered heterocyclie. cyano. or acts as 
the functional equivalent of I in the molecular recognition domam when bound to a 
5 TR, 

wherein R4' is 

-H, halogen, -CF3, -OH, -NH^, NH3, -NCCH^),, carboxylate, phosphonate. 
10 phosphate or sulfate, -SH, -CH3, -Et, or akyl, aryl or 5- or 6-membered 

heterocyclic aromatic attached through urea or carbamate linkages to O or N or S 
at the R,' position, or acts as the functional equivalent of OH in the molecular 
recognition domain when bound to a TR, 

15 wherein R5* is 

-H. -OH, -NH,, -N(CH3), -SH -NH3, -N(CH3)3. carboxylate. phosphonate, 
phosphate, sulfate, branched or straight chain alkyl having 1 to 9 carbons, 
substituted or unsubstituled aryl, wherein said substimted aryl is substituted with 

20 halogen or 1 to 5 carbon alkyl and wherein said aryl is optionally connected to the 

ring by a -CH^-. aromatic heterocycle having 5 to 6 atoms, wherein said 
heterocycle may be substituted with one or more groups selected from -OH, -NH,. 
-SH, -NHj, -N(CH3)3, carboxylate, phosphonate, phosphate or sulfate, heteroalkyl, 
arylalk^l. heteroaryl alkyl, polyaromatic, or polyheteroaromaiic. wherein said R,' 

25 may be substituted with polar or charged groups, 

wherein R^' is 

-H. halogen. -CF3, -OH, -NH,, -SH, -CH3, -Et, or acts as the functional 
30 equivalent of H in the molecular recognition domain when bound to a TR, 



wherein X is 
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O. S, so,. NH. NR;, CH,, CHR^. CR,R,, wherein R^ is alkyi, aryl or 5- or 
6-membered heterocyclic aromatic, 

and wherein said TR LBD hgand has an apparent Kd for binding TR LBD of 1 ^iM or 
5 less. 

65. The compound of claim 64, wherein 

Rj is carboxylate, phosphonate, phosphate or sulfite and is connected to the ring with a 
0 to 3 carbon linker, 
10 R2 IS H, 

R3 is -I, -Br, or -CHj, 
R5 IS -I, -Br, or -CH3. 
R« IS -H, 
R/ is -H, 

15 R3' is -I, -Br, -CH3, -iPr, -phenyl, benzyl, or 5- or 6-membered ring heterocyclcs, 
R,' is -OH, -NH2, and -SH, 

R5' is -H. -OH, -NHj, -N(CH3)2 -SH -NH^, -N(CH3)3, carboxylate, phosphonaie, 
phosphate, sulfate, branched or straight chain alkyl having 1 to 9 carbons, 
substituted or unsubstimted aryl, wherein said substituted aryl is substituted with 

20 halogen or 1 to 5 carbon alkyl and wherein said aryl is optionally connected to the 

ring by a -CHj-, aromatic heterocycle having 5 to 6 atoms, wherein said 
heterocycle may be substiaited with one or more groups selected from -OH, -NH;, 
-SH, -NH3. -NCCHj),, carboxylate, phosphonate, phosphate or sulfate, heteroalkyl, 
arylalkyl, heteroaryl alkyl, polyaromatic, or polyheteroaromatic, wherein said R5' 

25 may be substituted with polar or charged groups, and 

R^' is H. 

66 The compound of claim 65, wherem said compound is made by the method of 
claim 13. 

30 67. The compound of claim 65, wherein said compound is made by the method of 
claim 39. 

68. The compound of claim 65, wherein said compound is a TR antagonist. 
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69. The compound of claim 65. wherein said compound is a TR agomsi. 

70. The compound of claim 65, wherein said compound is a TR a selective Iigand. 

71. The compound of claim 65, wherein said compound is a TR ^ selective ligand. 

72. A pharmaceutical composition having selective thyromimetic activity comprising a 
5 compound as claimed in claim 65 and a pharmaceutically effective carrier. 

73. The composition of claim 72, wherein R5' is an alkyl having from 1 to 9 carbons 
and being straight chain or branched. 

74. A method of lowering the ratio of LDL-choiesterol to HDL-cholesterol levels 
comprising admimstering an effective amount of the compound of claim 65. 

10 75. A method of lowering plasma lipid levels which comprises administering to an 
animal in need thereof an effective amount of the compound of claim 65. 
76. A method of treating thyroid hormone deficiency in patients with compromised 
cardiac function, comprising administering an effective amount of the compound of claim 
65. 
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Figure 6 
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Figure 8 
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Dose-response to TS-10 -/+ 1 nM T3 in the 
TRAFalfal reporter cells 
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Dose-response to the compound TS-10 W- 1 nM T3 in the 
TRAFbetal reporter cell line 
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FIG. 19 





wo 97/21993 ^5/32 PCT^S96/20778 



Dose-response to TS-10 InM T3 in the liver 
reporter cell line HepG2 
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